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ABSTRACT 


Aspidoscelis tesselata exhibits significant clonal diversity despite its recent origin (from hybrid- 
ization between A. tigris marmorata and A. gularis septemvittata) and its parthenogenetic mode 
of reproduction. Two hypotheses have been advanced to explain the derivation of its genetic and 
morphological variation: (1) separate parthenogenetic lineages derived from several different F, 
hybrid zygotes, and (2) postformational mutations occurring in a parthenogenetic lineage derived 
from a single F, hybrid zygote. We evaluated these competing hypotheses with evidence from 
skin transplant studies, protein electrophoresis, multivariate analyses of morphological characters, 
and geographic distributions of pertinent groups. Starting with the clonal diversity at Conchas Lake 
State Park, San Miguel County, New Mexico, we expanded the study to include populations at 
Sumner Lake State Park and Fort Sumner (De Baca County), Puerto de Luna (Guadalupe County), 
and Arroyo del Macho and Roswell (Chaves County). This enabled us to resolve origins of color 
pattern classes and genotypic clones in eastern New Mexico. We used pattern class designations 
C-E and E-C to signify that elements of both pattern classes were expressed in populations at 
Conchas Lake and Arroyo del Macho. 

The two pattern classes at Conchas Lake (C-E and D) had the same F, hybrid karyotype 
(2n = 46), with haploid sets of 23 chromosomes characteristic of each progenitor species of 
A. tesselata. Clonal variation was found at 4 of the 35 gene loci examined electrophoretically: 
GPI (glucose-6-phosphate isomerase), EST2 (a muscle esterase), sACOH (aconitase hydra- 
tase), and MPI (mannose-6-phosphate isomerase). The strong congruence between genotype 
and morphological variation facilitated the characterization of three morphological subgroups 
of C-E. Although these subgroups lacked individually distinctive color patterns, they were 
discriminated effectively in canonical variate analyses based on scalation characters and a 
priori groups of known genotype. 

Nine individuals of Conchas C-E and four individuals of Conchas D have histocompatibility 
data from a recent skin transplant study (Cordes and Walker, 2003). The subgroup identities 
of the C-E specimens document histocompatibility among the three morphological subgroups 
of C-E and between each subgroup and representatives of pattern class D. This evidence, 
together with Maslin’s (1967) report of histocompatibility between pattern classes C and E, 
suggests that all color pattern classes, morphological subgroups, and genotypic clones of A. 
tesselata can be traced back to a single ancestral F, hybrid zygote. 

A pair of pale broken lines in the middorsal region distinguishes pattern class D from the other 
pattern classes. However, Conchas ID shared the GPI —100/—96, EST2 100/96 genotype with 
Conchas IC-E, and individuals of these pattern classes were very similar in multivariate meristic 
characters. Sumner D expressed the same type of relationship, resembling the syntopic population 
of Sumner C rather than the other population of D. In addition, certain individuals of Sumner C 
had partially divided (D-like) vertebral lines—additional evidence that Sumner C was ancestral to 
Sumner D. We conclude that pattern class New Mexico D is polyphyletic, having originated twice 
from different individuals of C-E and C in the vicinities of Conchas and Sumner Lakes. 

The northern position of pattern classes C and C-E in the range of A. tesselata is consistent 
with recent colonizations by individuals from more southerly populations. A candidate source 
population, based on its extensive color pattern and meristic variation, is E-C at Arroyo del 
Macho. The strong morphological resemblance of several northern populations to Macho E-C 
rather than to either syntopic clones or geographically proximate populations of other pattern 
classes supports this possibility. Evidence from geographic distributions, patterns of genotypic 
and meristic variation, and histocompatibility identifies postformational mutations as the likely 
basis for the genetic and morphological variation found in A. tesselata. This variation also 
includes different life-history characteristics between pattern classes C and E at Sumner Lake 
State Park. 

The name ftesselata is presently associated indirectly with pattern class C through the neo- 
type of A. tesselata. The neotype is a specimen of Colorado D, a derivative of pattern class 
C. With respect to pattern classes E-C, E, and other southern variants, taxonomic restructuring 
would confront mosaic patterns of genotypic, phenotypic, and geographic variation—patterns 
expected from random mutations in clonally reproducing species. Aspidoscelis tesselata has 
exploited a variety of ecological opportunities despite the constraints of clonal inheritance. 
Postformational mutations in the generalized genotype acquired from its progenitor species 
may have contributed to its ecological success. 
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INTRODUCTION 


Reeder et al. (2002) presented evidence 
that North American lizards previously as- 
signed to the genus Cnemidophorus repre- 
sent a monophyletic group, leaving South 
American species in possession of the ge- 
neric name in a separate clade. This neces- 
sitated resurrecting the generic name Aspi- 
doscelis for North American species. There- 
fore, the present paper involves an interpre- 
tation of the evolutionary history of the 
diploid parthenogenetic species Aspidoscelis 
tesselata in eastern New Mexico, the species 
formerly allocated to Cnemidophorus tesse- 
latus. 

Aspidoscelis tesselata is an excellent spe- 
cies for assessing the evolutionary potential 
of a parthenogenetic vertebrate. Its attributes 
include a broad latitudinal distribution, well- 
marked genotypic and phenotypic variation, 
and histocompatibility evidence that evolu- 
tionary divergence began with a single basal 
parthenogenetic individual (Maslin, 1967; 
Cordes and Walker, 2003). Of additional evo- 
lutionary interest, Aspidoscelis tesselata par- 
ticipated in the origin of a new species. Al- 
though hybridization between A. tesselata 
and A. figris marmorata produced sterile 
triploid hybrids (Taylor et al., 2001), the trip- 
loid species A. neotesselata (Walker et al., 
1997a) originated from an A. tesselata X A. 
sexlineata hybridization (Parker and Selan- 
der, 1976; Densmore et al., 1989; Dessauer 
and Cole, 1989). Equally intriguing is the 
possibility that a new parthenogenetic spe- 
cies could result from genetic changes (di- 
vergence) in an established parthenogenetic 
species (Echelle, 1990; Taylor and Cooley, 
1995a, 1995b; Manriquez-Mordn et al., 
2000; Schmitz et al., 2001). In addition to 
this possibility, we investigated the origin of 
color pattern classes and the pattern of di- 
vergence that has occurred in A. tesselata 
from eastern New Mexico. 

Aspidoscelis tesselata originated recently 
(Densmore et al., 1989; Reeder et al., 2002) 
from hybridization between a female A. tigris 
marmorata and a male A. gularis septemvit- 
tata (Neaves, 1969; Parker and Selander, 
1976; Dessauer and Cole, 1989; Dessauer et 
al., 1996). Despite its clonal pattern of in- 
heritance (Dessauer and Cole, 1986), A. tes- 
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selata is ecologically successful based on a 
geographic range that includes parts of Chi- 
huahua, Mexico (Zweifel, 1965), Texas (Dix- 
on, 2000), New Mexico (Degenhardt et al., 
1996), Oklahoma (Webb, 1970), and Colo- 
rado (Hammerson, 1999), and its presence in 
various assemblages of bisexual and parthe- 
nogenetic congeners (Cuellar, 1979; Taylor et 
al., 2000, 2001). Based on its extensive lat- 
itudinal distribution, Aspidoscelis tesselata 
might have inherited an advantageous gen- 
eralized genotype from its two progenitor 
species (Taylor et al., 2001). This ecological 
success has occurred despite the absence of 
genetic recombination. Price (1992) and 
Price et al. (1993) provide a contrasting in- 
terpretation of “‘success”’ in A. tesselata. 

Morphological variation in A. tesselata in- 
cludes four color pattern classes (Zweifel, 
1965; Taylor et al., 1996) integrated into the 
currently understood phylogenetic history of 
A. tesselata by Walker et al. (1997a). Anoth- 
er diploid pattern class (F), included in A. 
tesselata by Zweifel (1965), was described 
as A. dixoni by Scudday (1973), and Dens- 
more et al. (1989) presented mitochondrial 
DNA evidence that A. dixoni and A. tesselata 
may have originated from independent hy- 
bridizations. 

Color pattern differences among pattern 
classes C, Colorado D, New Mexico D, and 
E were used to discover local examples of 
ecological (Walker et al., 1997b) and life-his- 
tory (Taylor et al., 1997, 2000) differences 
between sympatric clones. However, we ex- 
pose (below) complications in using color 
pattern classes C and E as simple descriptors 
of regional patterns of geographic variation. 

Our study is rooted in the important pio- 
neering investigations of Zweifel (1965), 
Maslin (1967), Parker and Selander (1976), 
and Parker (1979). Zweifel (1965) set the 
stage for future research by recognizing that 
a few color pattern classes could accommo- 
date the previously undecipherable morpho- 
logical variation in A. tesselata. He was then 
able to describe A. tesselata from Conchas 
Lake as comprising two color pattern classes 
(C and D) and document the remarkable me- 
ristic variation in Conchas C. The magnitude 
of this morphological variation was unex- 
pected because it challenged the intuitive ex- 
pectation that parthenogenetic reproduction 
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constrains phenotypic variability. With great 
acuity, Zweifel (1965) predicted that Con- 
chas C might consist of several clones, and 
electrophoretic analyses of proteins by Park- 
er and Selander (1976) confirmed this pre- 
diction. Parker (1979) then demonstrated 
congruence between meristic and genotypic 
variation, thereby providing genetic markers 
(clones IC, VIC, and VIIIC) for morpholog- 
ical subgroups of Conchas C. These clones 
appeared to have combinations of alleles that 
were polymorphic in their bisexual ancestors. 
Therefore, Parker and Selander (1976), Park- 
er (1979), and Parker et al. (1989) argued 
that several independent hybridizations were 
responsible for the genetic and morphologi- 
cal variation in Conchas C. However, Maslin 
(1967) revealed histocompatibility between 
pattern classes C and E, implying that the 
genealogies of different pattern classes con- 
verged on a single parthenogenetic progeni- 
tor. 

The purpose of the present study was to 
(1) reconcile patterns of genotypic and mor- 
phological variation for clones of A. tesselata 
at Conchas Lake and different color pattern 
classes from eastern New Mexico, (2) reas- 
sess the multiple hybridization hypothesis to 
explain variation among these groups, (3) 
evaluate evidence that the genetic and mor- 
phological variation in A. tesselata was gen- 
erated by postformational mutation, and (4) 
determine if divergent clones deserve formal 
taxonomic recognition. 


MATERIALS AND METHODS 
MORPHOLOGICAL ANALYSES 


In a departure from all previous studies, 
we use new pattern class designations to flag 
color pattern ambiguities in A. tesselata at 
Conchas Lake State Park and Arroyo del Ma- 
cho. The results of a discriminant analysis 
(DA) justified two pattern class designations. 
Using Zweifel’s (1965) pattern class desig- 
nations, the DA model was based on samples 
of pattern class “‘C”’ from Conchas Lake and 
pattern class “‘E”’ from Arroyo del Macho as 
the a priori groups (fig. 1). These specimens 
met body size criteria (see below) that re- 
moved ontogenetic variation from the three 
color pattern characters used in the DA (table 
1). The DA model, with standardized dis- 
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criminant functions of 1.308 for L-breaks, 
—0.519 for DL-breaks, and —0.380 for PV- 
breaks, classified 26 of 100 specimens of 
Conchas “‘C”? as Macho “E”’ and 14 of 37 
specimens of Macho “E”’ as Conchas “‘C’’. 
Despite misclassified individuals and a weak 
eigenvalue (0.113), the two groups differed 
significantly in discriminant function scores 
(t135 = —3.912; P < 0.0005). Therefore, we 
used two new pattern class designations: C-E 
for the Conchas Lake non-Ds, and E-C for 
the form at Arroyo del Macho. The first letter 
signifies the pattern class identity of this pop- 
ulation in previous studies; the second letter 
emphasizes that this population contains a 
range of color pattern variation that includes 
individuals with color pattern features of this 
second pattern class. Additional evidence of 
color pattern ambiguities in these two groups 
is presented in other sections of the paper. 
With the addition of two new color pattern 
designations, our study was based on five 
color pattern classes of A. tesselata from six 
localities. The focal samples comprised pat- 
tern classes C-E (N = 135) and D (N = 51) 
from Conchas Lake State Park (fig. 1). These 
samples included 34 specimens (21 C-E, 13 
D) with protein phenotypes identified elec- 
trophoretically by Parker and Selander 
(1976) and H.C.D. and C.J.C. (this study) 
and 13 specimens (9 C-E, 4 D) used in skin 
transplant experiments by Cordes and Walker 
(2003). Two additional groups of samples 
were included to determine patterns of evo- 
lutionary divergence. The first group con- 
sisted of pattern classes E and E-C with elec- 
trophoretic data from Fort Sumner (N = 11), 
Roswell (NV = 10), and Arroyo del Macho (NV 
= 10). The second group comprised samples 
lacking electrophoretic data. These were (1) 
pattern class C (N = 47), pattern class E (NV 
= 47), and pattern class D (N = 5) from 
Sumner Lake State Park (approximately 90 
linear km south of the Conchas Lake locali- 
ties), (2) pattern class E (N = 25) from 
Puerto de Luna (approximately 34 linear km 
northwest of the Sumner Lake collecting 
sites), (3) pattern class E (N = 5) from Fort 
Sumner (approximately 18 linear km south- 
east of Sumner Lake), and (4) pattern class 
E-C (N = 31) from Arroyo del Macho (ap- 
proximately 105 linear km south of Sumner 
Lake). The geographic relationships of sam- 
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Fig. 1. Geographic relationships among four northern collecting localities of Aspidoscelis tesselata 
of color pattern classes C, New Mexico D, and E and convenience classes C-E and E-C. Color patterns 
found at the four sites are (1) Conchas Lake State Park: C-E and New Mexico D; (2) Sumner Lake 
State Park: C, New Mexico D, and E; (3) Puerto de Luna: E; and (4) Arroyo del Macho: E-C. 
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TABLE 1 
Meristic Characters (and SVL) Used in Analyses of Morphological Variation in 
Aspidoscelis tesselata from the Locations Studied 














Character 
abbreviation Description 

GAB Number of granular dorsal scales in a single row around midbody. There are eight longitudinal rows of 
enlarged scales making up the ventral body surface. The third ventral row on either side of the mid- 
sagittal line terminates anteriorly in the axillary region. The 15th ventral scale posterior to this terminus 
established the point for beginning the GAB count. 

COS Bilateral total of circumorbital scales as standardized by Wright and Lowe (1967). 

LSG Sum of lateral supraocular granules on both sides of the head. These granular scales are located between 
the supraoculars and superciliary scales, and the count includes all scales anterior to the suture line 
between the third and fourth supraoculars. 

FP Sum of femoral pores on both thighs. 

GS Number of granular gular scales bordering the medial edges of the eight anterior sublabials (four on each 
side, or their subdivisions) and the posterior mental. 

PSC Sum of all scales, including occipitals, contacting the outer perimeter of parietal and interparietal scales. 

L-breaks Total number of interruptions by black pigment of the lateral pale stripes. 

DL-breaks Total number of interruptions by black pigment of the dorsolateral pale stripes. 

PV-breaks Total number of interruptions by black pigment of the paravertebral pale stripes. 

SDL-T4 Number of subdigital lamellae on the fourth toe of one foot (right was chosen unless damaged). 

SVL Length of body from tip of snout to posterior edge of preanal scales (in mm). 


pling localities are shown in figure 1, and 
specific information on samples and sam- 
pling localities is provided in appendix 1. 

We shortened frequent references to geno- 
typic clones (GPI and EST2 loci) and affili- 
ated pattern class as follows: Conchas IC-E 
(Roman numerals indicate that the genotype 
is known) and 1C-E (Arabic numerals refer 
to morphological subgroups corresponding to 
a priori groups of known genotype); VIC-E 
and 6C-E; VHIC-E and 8C-E. We also com- 
bined abbreviated names of collecting local- 
ities and pattern class designations to refer to 
various groups (Conchas D; Sumner C, D, 
and E; Luna E; Macho E-C). 

Each specimen was scored for 10 meristic 
characters representing seven scalation fea- 
tures and three color pattern features (tables 
1, 2), the latter involving the number of 
breaks (disruptions by transverse black bars) 
in the lateral, dorsolateral, and paravertebral 
pale stripes. A traditional character (SDL-T4 
in tables 1 and 2) was omitted from multi- 
variate analyses to maintain consistency and 
comparability across different color pattern 
classes, genotypic clones, and morphological 
subgroups. We excluded this character be- 





cause certain individuals in critical samples 
(those with protein phenotypes and those 
used for skin transplants) had damaged or 
missing 4th toes on both feet. 

Body length (snout—vent length = SVL) 
was measured to the nearest millimeter with 
digital calipers. Significant relationships 
were found between color pattern characters 
and SVL in composite samples of pattern 
classes C and C-E (N = 173) and E and E-C 
(N = 112). These relationships were between 
L-breaks and SVL in C + C-E (P = 0.003) 
and between L-breaks and SVL (P = 0.016) 
and PV-breaks and SVL (P = 0.001) in E + 
E-C. Using specimens >66 mm SVL for pat- 
tern classes C, C-E, and D and >69 mm SVL 
for pattern classes E and E-C removed on- 
togenetic variation from the definitive sam- 
ples (P = 0.08 and P = 0.62 for the regres- 
sions of L-breaks on SVL, and P = 0.18 for 
the regression of PV-breaks on SVL). 

Three morphological subgroups of pattern 
class C-E were identified from a canonical 
variate analysis (CVA) of 21 individuals of 
known genotype, using the three GPI, EST2 
genotypes as a priori groups. These sub- 
groups, corresponding to genotypic clones, 
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were then used as a priori groups to assign 
individuals of C-E lacking electrophoretic 
data to morphological subgroup. The pattern 
of meristic variation among the morpholog- 
ical subgroups was then depicted by a defin- 
itive CVA, with subgroup assignments used 
to define a priori groups. The same procedure 
was used to assign the skin transplant spec- 
imens of Cordes and Walker (2003) to mor- 
phological subgroup and show their patterns 
of variation. 

Nine meristic characters (table 1) were 
used in most analyses. Exceptions included 
CVAs based only on color pattern characters 
and those used for skin transplant specimens. 
We excluded characters based on DL-breaks 
and PV-breaks from the latter because the 
dorsolateral and paravertebral stripes had 
been used as transplantation sites, and the in- 
tegrity of the pattern was compromised in 
these specimens. Only specimens with com- 
plete data were included in each analysis. 
Specific procedures are detailed in the appro- 
priate sections of Results and Discussion. 

To facilitate comparisons, a symmetric dis- 
similarity matrix of Mahalanobis distances 
(D’) between group means was used to con- 
struct an additive tree depicting the morpho- 
logical similarities and differences among the 
various groups. The D? values were obtained 
from BMDP program 5M, Linear and Qua- 
dratic Discriminant Analysis, accessed via 
SYSTAT 10.2. Additive trees are directed 
graphs with paths (sums of horizontal branch 
lengths) representing relative distances be- 
tween entities. This procedure, in conjunc- 
tion with ordination (de Queiroz and Good, 
1997), is more appropriate than hierarchical 
clustering for illustrating patterns of diver- 
gence from similarity/dissimilarity data (Sat- 
tath and Tversky, 1977; Wilkinson et al., 
1996; de Queiroz, 1998). For example, intra- 
cluster branches (representing specific 
groups) can vary in length on an additive tree 
but not on trees generated by hierarchical 
joining techniques. Therefore, additive trees 
provide information on differences among 
intracluster groups in addition to intercluster 
groups. Additive trees are not intended as 
phylogenetic reconstructions, but clonal di- 
vergence manifested in quantitative meristic 
differences should be evident. 


TAYLOR ET AL.: ASPIDOSCELIS TESSELATA 9 


ELECTROPHORETIC ANALYSES 


Two of the 21 presumptive protein loci as- 
sayed with horizontal starch gel electropho- 
resis by Parker and Selander (1976), Pgi (= 
GPI, glucose-6-phosphate isomerase) and 
Est-2 (= EST2, a muscle esterase), and 4 of 
the 31 presumptive protein loci assayed by 
H.C.D. and C.J.C. (GPI, EST2, aconitase hy- 
dratase [SACOH], and mannose-6-phosphate 
isomerase [MPI]) revealed clonal differences 
in the Conchas Lake samples. Parker and Se- 
lander (1976) designated alleles at each locus 
according to migration of their allozymes rel- 
ative to the most common allele (“‘100”’ for 
anodally migrating proteins and “‘—100” for 
cathodally migrating proteins) found in A. ti- 
gris marmorata. Clonal designations from 
Parker (1979) use a combination of multilo- 
cus electrophoretic genotype (Roman numer- 
als I-XIII) and pattern class (B = A. neotes- 
selata; C, D, and E = A tesselata; F = A. 
dixoni). All available specimens of pattern 
classes C (= C-E here) and D analyzed by 
Parker and Selander (1976) and Parker 
(1979) were rescored for meristic characters 
and combined with AMNH specimens with 
known GPI and EST2 genotypes from Con- 
chas Lake to represent the foundation for the 
present study: 10 individuals of clone IC-E 
(GPI —100/—96, EST2 100/96), 7 individu- 
als of clone VIC-E (GPI —99/—96, EST2 
100/96), 4 individuals of clone VUIC-E (GPI 
—100/—96, EST2 96/96), and 10 individuals 
of clone ID (GPI —100/—96, EST2 100/96). 
Unfortunately, the three individuals of clone 
VIIC (GPI —99/—96, EST2 96/96) described 
by Parker (1979) have been lost and were 
thus unavailable for reanalysis. 

New electrophoretic data for 31 gene loci 
have been provided by H.C.D. and C.J.C. 
(table 3), based on more recently collected 
AMNH specimens (see appendix 1). These 
loci include 17 of the 21 analyzed by Parker 
and Selander (1976). Methodology followed 
Harris and Hopkinson (1976), Murphy et al. 
(1996), and, particularly for North American 
lizards of the genus Aspidoscelis, Dessauer 
et al. (2000). For each locus, alleles are des- 
ignated in alphabetical order according to de- 
creasing anodal migration of their allozymes. 
For multilocus enzymes, loci are listed nu- 


10 AMERICAN MUSEUM NOVITATES NO. 3424 


TABLE 3 
Genotypes at 31 Gene Loci* in Clones? of Aspidoscelis tesselata of Pattern Classes C—E and D 
from Conchas Lake State Park Compared with Pattern Class E-—C from Arroyo del Macho 
(See Taylor et al., 2001.) 








TESC-E TESC-E TESC-E TESD TESE-C 
Locus (N = 3) (N = 2) (N = 2) (N = 2) (N = 11) 
Oxidoreductases 
ADH ab ab ab ab ab 
G3PDH aa aa aa aa aa 
IDDH aa aa aa aa aa 
LDHI1 ab ab ab ab ab 
LDH2 aa aa aa aa aa 
sMDH ab ab ab ab ab 
mMDH aa aa aa aa aa 
sMDHP ab ab ab ab ab 
sIDH bb bb bb bb bb 
mIDH aa aa aa aa aa 
sSOD ab ab ab ab ab 
mSOD aa aa aa aa aa 
Transferases 
sAAT ab ab ab ab ab 
mAAT ab ab ab ab ab 
CK1 aa aa aa aa aa 
AK aa aa aa aa aa 
Hydrolases 
ESTD bb bb bb bb bb 
EST2¢ be be bc be be 
PEPA cd cd cd cd cd 
PEPB bb bb bb bb bb 
PEPD cc cc cc cc cc 
PEPE aa aa aa aa aa 
ADA ac ac ac ac ac 
Lyase 
sACOH ac be be be be 
Isomerases 
MPI ac ac ab ac ab 
GPlIe ac ac ab ac ac 
PGM2 ab ab ab ab ab 
PGM3 ab ab ab ab ab 
Blood proteins 
TFe bb bb bb bb bb 
ALB ab ab ab ab ab 
HB aa aa aa aa aa 








4Abbreviations for loci are as follows: ADH, alcohol dehydrogenase; G3PDH, glycerol-3-phosphate dehydrogenase; IDDH, 
L-iditol dehydrogenase; LDH, L-lactate dehydrogenase; MDH, malate dehydrogenase; MDHP, malate enzyme; IDH, isocitrate dehy- 
drogenase; SOD, superoxide dismutase; AAT, aspartate aminotransferase; CK, creatine kinase; AK, adenylate kinase; EST, esterase; 
PEP, peptidase; ADA, adenosine deaminase; ACOH, aconitase hydratase; MPI, mannose-6-phosphate isomerase; GPI, glucose-6- 
phosphate isomerase; PGM, phosphoglucomutase; TF, transferrin; ALB, albumin; and HB, hemoglobin; s, cytosolic enzyme; m, 
mitochondrial enzyme. For loci at which no a- or b-allele is shown, these occurred in other samples or in closely related taxa that 
are not included in this paper. 

bDistinctive alleles (not found in pattern class E-C) that characterize different clones of pattern classes C-E and D are in an 
enlarged, boldface type and involve only three loci: sACOH, MPI, and GPI. 

©This table consists of new data determined by H.C.D. and C.J.C. For EST2, all lizards were heterozygotes, thus consistent with 
the common clone designated by Parker and Selander (1976) and Parker (1979) as 100/96. For GPI, the ac genotype is the same 
as the common clone designated as — 100/—96 for Pgi by Parker and Selander (1976) and Parker (1979), whereas our ab geno- 
type is their —99/—96 genotype. As we only recently recognized the subtle differences in GPI, the ac genotype we report here is 
the same as the ab genotype reported by Taylor et al. (2001). 
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merically in order of decreasing anodal mi- 
gration of their isozymes. 


KARYOTYPIC ANALYSIS 


We used previously published methods for 
preparing and studying standard, giemsa- 
stained chromosomes (Cole, 1979). We ex- 
amined 38 cells at mitotic metaphase from 
bone marrow of four females of pattern class 
C-E and three individuals of pattern class D 
of A. tesselata from Conchas Lake State 
Park, San Miguel County, New Mexico (see 
appendix 1). 


RESULTS AND DISCUSSION 


COLOR PATTERN FEATURES OF PATTERN 
CLASSES ARE INHERITED CLONALLY 


Partitioning of morphological variation in 
A. tesselata traditionally begins with an as- 
signment of individuals to color pattern class. 
Therefore, it is important to know that fun- 
damental elements of color pattern are trans- 
ferred clonally from mothers to daughters. 
Dessauer and Cole (1986) established this 
fact for pattern classes C-E and D from Con- 
chas Lake State Park, but not all individuals 
documenting this clonal inheritance were 
specified. The evidence included three adults 
of pattern class C-E from Conchas Lake that 
produced eggs in the Aspidoscelis colony at 
the American Museum of Natural History. 
The nine hatchlings derived from these 
clutches were all non-D in pattern class; that 
is, they lacked the paired, broken vertebral 
lines that characterize pattern class D 
throughout ontogeny. Similarly, two repre- 
sentatives of pattern class D from Conchas 
Lake produced eight hatchlings from two 
clutches of eggs; all had the characteristic 
color pattern of New Mexico D. All sets of 
hatchlings and their mothers are identified in 
appendix 1. 


KARYOTYPES ARE CONSISTENT WITH 
HYBRID ORIGIN 


To date, two clearly resolved diploid kar- 
yotypes have been published for Aspidoscelis 
tesselata. Dessauer and Cole (1989: 57, fig. 
8A) depict one from an individual of pattern 
class D from Conchas Lake State Park, and 
Taylor et al. (2001: 24, fig. 1OA) depict one 
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from an individual of pattern class E-C from 
Arroyo del Macho. These two karyotypes are 
identical, and as discussed in detail by Taylor 
et al. (2001), the karyotypes illustrated are 
identical to the expectations for F, hybrids 
from A. tigris marmorata X A. gularis sep- 
temvittata (a haploid set of chromosomes 
from each; 2n = 46). A variant karyotypic 
clone was found at Arroyo del Macho (also 
pattern class E-C) in which the X-chromo- 
some originally inherited from A. tigris mar- 
morata was apparently fissioned at the cen- 
tromere to become two smaller telocentric 
chromosomes (2n = 47). 

For the present study, we determined the 
karyotypes of four specimens of A. tesselata 
pattern class C-E and three specimens of pat- 
tern class D from Conchas Lake State Park. 
All individuals had the typical F, hybrid kar- 
yotype (2n = 46) with the original, unfis- 
sioned X-chromosome from A. figris. 


GENOTYPIC DIFFERENCES WITHIN PATTERN 
CLASS C-E FROM CONCHAS LAKE 


PARKER AND SELANDER ANALYSES: Allo- 
zyme variation at 21 gene loci was reported 
in a range-wide study of A. tesselata (includ- 
ing the triploid A. neotesselata) by Parker 
and Selander (1976), who identified process- 
es that could account for the patterns of ge- 
notypic variation detected (Parker, 1979; 
Parker et al., 1989). The three representatives 
of clone VIIC have been lost and cannot be 
reanalyzed morphologically for canonical 
variate analysis. Clone VIIC shares the —99/ 
—96 GPI genotype with clone VIC and the 
96/96 EST2 genotype with clone VIIIC 
(Parker, 1979). One individual of clone VIIC 
was similar to clones IC-E, VHIC-E, and ID 
for the nine meristic characters scored by 
Parker (1979), while the other two individ- 
uals were similar to clone VIC-E in having 
significantly higher counts for LCOF (eft 
circumorbital scales in contact with frontal 
scale), GAB (granules around midbody at the 
level of the 20th ventral scale row, as scored 
by Parker, 1979), L3 (number of lamellae on 
the 3rd toe of the left hindfoot), and SAT 
(number of scales around the tail at the level 
of the 15th caudal scale row). Thus, eight of 
the nine individuals of clones VIC-E and 
VIUC-E, marked by the GPI —99/—96 ge- 
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TESE (ac) 
TESC (ab) 
TESC (ac) 
TESC (ab) 
TESC (ac) 
TESC (ac) 





cba 


Fig. 2. Electrophoretic phenotypes of GPI, a 
dimeric enzyme, from erythrocyte hemolysates of 
six specimens of Aspidoscelis. Letters below gel 
identify allozymes based on alleles present (table 
3), and the genotype of each lizard is listed on 
the right. Note the very slight difference in mi- 
gration between the products of the b-allele versus 
c-allele. Lanes for individual lizards are labeled 
beside their patterns on the gel as follows: TESC, 
A. tesselata of pattern class C-E from Conchas 
Lake State Park, New Mexico; and TESE, A. tes- 
selata of pattern class E from Sandoval County, 
New Mexico. Anode is to the right. 


notype (Parker, 1979, table 2: 1153), were 
morphologically divergent from clones 
marked by GPI —100/—96. Although clone 
VIIC is apparently rare, some of its repre- 
sentatives may have been included in our 
morphological subgroups. However, we are 
confident that the patterns of variation have 
not been distorted because clone VIIC is not 
distinctive in meristic variation (Parker, 
1979). 

DESSAUER AND COLE ANALYSES: Seven 
specimens of pattern class C-E and two of 
pattern class D from Conchas Lake State 
Park were assessed at 31 gene loci (table 3). 
As was found previously (Parker and Selan- 
der, 1976), GPI genotypes of —100/—96 (IC- 
E) and —99/—96 (VIC-E) were present (ac 
and ab for GPI, respectively, in table 3; fig. 
2), but all nine individuals were apparently 
heterozygous for the EST2 alleles (bc for 
EST2 in table 3, presumably the same as 
Parker and Selander’s 100/96). Note, how- 
ever, that the gene products revealed at EST2 
may not actually represent only one gene lo- 
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TESC (ac) 
TESC (ac) 
TESC (bc) 
TESC (bc) 
TESC (bc) 
TESC (ac) 
TESC (bc) 
TESC (ac) 
TESC (be) 
c ba 
Fig. 3. Electrophoretic phenotypes of 


sACOH, a monomeric enzyme, from liver ho- 
mogenates of nine specimens of A. tesselata of 
pattern class C-E from Conchas Lake State Park, 
New Mexico. Letters below gel identify allo- 
zymes based on alleles present (table 3), and the 
genotype of each lizard is listed on the right. 
Lanes for individual lizards are labeled beside 
their patterns on the gel. Anode is to the right. 


cus, as they do with reasonable certainty for 
the other loci tested. In addition, genotypic 
variation was detected at two other loci (table 
3; sACOH, fig. 3, and MPI, fig. 4), neither 
of which was included in the Parker and Se- 
lander (1976) analyses. This additional ge- 
notypic variation included genotypes be and 
ac at the sACOH locus and genotypes ac and 
ab at the MPI locus. These genotypic des- 
ignations denote the relative migration dis- 
tances in electrophoretic gels, with the al- 
phabetical sequence of letters representing 
decreasing anodal migration. 

A review of all of the electrophoretic data 
presented for A. tesselata of pattern classes 
C, C-E, E-C, E, and New Mexico D by Park- 
er and Selander (1976), Dessauer and Cole 
(1986, 1989), Taylor et al. (2001), table 3, 
and additional unpublished data of H.C.D. 
and C.J.C. reveals that data are now available 
for a total of 34 clearly identified gene loci 
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NEOTESA _ (bbc) 
NEOTESB (abc) 
| NEOTESC (acc) 
TESC (ac) 
TESD (ac) 
TESE (ab) 
TESF (aa) 
TESF (aa) 


TESF x PUN (aaa) 
TESG (ab) 
| TESH (ab) 





c b a 


Fig. 4. Electrophoretic phenotypes of MPI, a 
monomeric enzyme, from liver homogenates of 
11 specimens of Aspidoscelis. Letters below gel 
identify allozymes based on alleles present (table 
3), and the genotype of each lizard is listed on 
the right. Lanes for individual lizards are labeled 
beside their patterns on the gel (with genotype) as 
follows: NEOTESA, B, and C, different pattern 
classes of the triploid A. neotesselata from Colo- 
rado; TESC and D, A. tesselata of pattern classes 
C-E and D from Conchas Lake State Park, New 
Mexico; TESE, A. tesselata of pattern class E-C 
from Arroyo del Macho, New Mexico; TESE A. 
dixoni from New Mexico; TESF xX PUN, triploid 
hybrid of A. dixoni X A. tigris punctilinealis from 
New Mexico; and TESG and H, A. dixoni of two 
pattern classes from Texas. Anode is to the right. 


(excluding the questionable EST2) for which 
comparisons have also been made with the 
bisexual parental taxa, A. tigris marmorata 
and A. gularis septemvittata. Heterozygosity 
is fixed and very high, at 47% of these loci, 
in the parthenogenetic A. tesselata, whereas 
it averages 5% for bisexual species of Aspi- 
doscelis (Parker and Selander, 1976; Des- 
sauer and Cole, 1989). As initially pointed 
out by Parker and Selander (1976), the gen- 
erality remains that for all loci the combi- 
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nation of alleles found in A. ftesselata is a 
combination that could be found in an F, hy- 
brid derived from A. tigris marmorata X A. 
gularis septemvittata, with very few excep- 
tions in variant clones (e.g., sACOH and 
MPI). Furthermore, the allele combinations 
present in individuals of pattern classes C, C- 
E, New Mexico D, E, and E-C typically are 
identical or nearly identical to each other for 
each locus. 


CONGRUENCE BETWEEN EXTERNAL 
MORPHOLOGY AND GENETIC MARKERS IN 
PATTERN CLASS C-E FROM CONCHAS LAKE 
STATE PARK 


Three subgroups of pattern class C-E at 
Conchas Lake are morphologically congru- 
ent with GPI, EST2, sACOH, and MPI ge- 
notypes (fig. 5). These subgroups were char- 
acterized by a canonical variate analysis of 
meristic characters, using the 21 specimens 
with known GPI and EST2 genotypes (Park- 
er’s [1979] clonal identities) as three a priori 
groups: clone IC-E (GPI —100/—96, EST2 
100/96), clone VIIIC-E (GPI -—100/—96, 
EST2 96/96), and clone VIC-E (GPI —99/ 
—96, EST2 100/96). We included the three 
color pattern characters in the CVA after first 
regressing each character on SVL and dem- 
onstrating an absence of ontogenetic varia- 
tion in the a priori groups (L-breaks: 7? = 
0.10, P = 0.17; DL-breaks: 7? < 0.0005, P 
= 0.98; PV-breaks: r? = 0.07, P = 0.24). The 
CVA model (table 4) had an overall classi- 
fication success of 90%, with the misclassi- 
fications being 2 of 10 individuals of clone 
IC-E classified to clone VIIIC-E (which we 
identify below as a likely mutational deriv- 
ative of clone IC-E). 


PATTERNS OF VARIATION AMONG 
MORPHOLOGICAL SUBGROUPS OF PATTERN 
CLASS C-E AND PATTERN CLASS D FROM 

CONCHAS LAKE STATE PARK 


Samples of known genotype (clones IC-E, 
N = 10; VIC-E, N = 7; VUIC-E, N = 4) 
were used as a priori groups in a CVA to 
assign to morphological subgroup those 
specimens of Conchas C-E with unknown 
genotypes but having complete meristic data 
(N = 79). Arabic numerals designated sub- 
groups (1C-E, 6C-E, and 8C-E) correspond- 
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ing to genotypic clones (IC-E, VIC-E, and 
VIIC-E). The CVA assigned 48 individuals 
to 1C-E, 18 to 6C-E, and 13 to 8C-E. The 
assignment CVA was followed by a second 
CVA (table 5), using subgroups 1C-E, 6C-E, 
and 8C-E as the a priori groups. This CVA 
had a classification success of 96%, which 
was reflected in the clear delineation of sub- 
groups in the ordination of canonical variate 
scores (fig. 6A). Three misclassifications 
were made among the 79 specimens—one 
individual of 6C-E assigned to 8C-E and two 
individuals of 1C-E, one assigned to each of 
subgroups 6C-E and 8C-E. Another CVA 
(table 6) substantiated that none of the sub- 
groups (1C-E, 6C-E, and 8C-E) had a dis- 
tinctive color pattern (fig. 6B). The classifi- 
cation success was only 42%, and each mor- 
phological subgroup contained individuals 
with C-like and E-like color pattern features 
(figs. 7, 13). 

To determine the morphological subgroup 
most closely resembled by Conchas D, a 
CVA (table 7) was based on a priori groups 
composed of morphological subgroups 1C-E, 
6C-E, 8C-E, and pattern class D. Although 
Parker and Selander (1976) reported single 
representatives of clones VID and VIID 
from Conchas Lake, the voucher specimen 
for clone VIIID (EDP 842) was reidentified 
aS a representative of pattern class C by 
Walker et al. (1997a) and clone VIIIC-E (this 
study). The voucher specimen for clone VID 
(EDP 846) is missing. However, Parker 
(1979) determined that it was morphologi- 
cally similar to individuals of clone ID. Our 
sample of ID comprises 21 specimens, in- 
cluding the 9 specimens analyzed electropho- 
retically by Parker and Selander (1976) and 
Parker (1979) and 1 specimen (EDP 835) 
identified as a representative of clone IC in 
those publications. Walker et al. (1997a) as- 
signed the latter specimen to pattern class D, 


< 
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TABLE 4 
Discriminant Functions Used to Distinguish 
Genotypically Different Groups of 
Aspidoscelis tesselata of Pattern Class C—E 
from the Vicinity of Conchas Lake, San 
Miguel County, New Mexico 











Character CV1 CV2 
GAB —1.215 0.202 
L-breaks 0.826 0.284 
LSG —0.340 0.008 
DL-breaks 0.089 ~0.975 
PV-breaks —0.321 0.668 
Cos 0.471 —0.522 
FP —0.213 —-0.411 
PSC 0.046 0.299 
GS 0.043 —0.265 
Eigenvalue 15.626 0.915 


Total explained variation 94.5% 5.5% 





«Characters are defined in table 1. 


and we concur, but these differences of opin- 
ion are expected to arise when visually clas- 
sifying specimens with ambiguous color pat- 
terns. Eleven specimens of Conchas D ana- 
lyzed electrophoretically by H.C.D. and 
C.J.C. were identical to clone IC-E for all 
loci tested and thus represented clone ID. Al- 
though it appears that clonal variation in pat- 
tern class D is limited, we use the more gen- 
eral designation of D for this a priori group. 

Color pattern differences between Con- 
chas C-E and Conchas D are centered in the 
vertebral field, where a pair of fragmented 
vertebral lines distinguishes the latter (figs. 
7, 8). This distinctive color pattern difference 
belies the fact that clones ID and IC-E are 
identical to each other in electrophoretic 
characters, and Conchas D and morphologi- 
cal subgroupIC-E are nearly identical in 
multivariate morphological characters (fig. 
9A, B). 


A. Pattern of multivariate morphological variation among 21 individuals of Aspidoscelis 


tesselata of pattern class C-E with genotypic data from Conchas Lake State Park, San Miguel County, 
New Mexico. Fourteen of these specimens were used by Parker and Selander (1976) to identify, from 
GPI and EST2 genotypes, clones IC-E, VIC-E, and VIIC-E. B. Open symbols denote Parker and 
Selander specimens shown in part A and solid symbols identify the seven AMNH specimens in part A 
with known GPI, EST2, sACOH (given first), and MPI (given second) genotypes. Canonical variate 
scores were derived from a canonical variate analysis (table 4) using nine meristic characters. 
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TABLE 5 
Discriminant Functions Used to Distinguish 
Three Morphological Subgroups of Pattern 
Class C-E of Aspidoscelis tesselata from the 
Vicinity of Conchas Lake, San Miguel 
County, New Mexico 


Character@ CV1 CV2 
GAB —0.959 —0.067 
COS 0.537 -0.311 
PV-breaks 0.118 0.052 
GS 0.111 0.097 
L-breaks 0.079 0.059 
FP 0.310 —0.790 
PSC 0.025 0.477 
DL-breaks 0.186 —0.310 
LSG —0.212 —0.186 
Eigenvalue 4.284 1.188 


Total explained variation 78.3% 21.7% 





4Characters are defined in table 1. 


HISTOCOMPATIBILITY AMONG 
MORPHOLOGICAL SUBGROUPS OF PATTERN 
CLASS C-E AND PATTERN CLASS D 


Cordes and Walker (2003) documented 
histocompatibility among the nine individu- 
als of Conchas C-E tested. We determined 
the morphological subgroup membership of 
each specimen by using CVA and individuals 
of known genotype (10 representatives of IC- 
E, 7 of VIC-E, and 4 of VIIIC-E) as a priori 
groups. The CVA model assigned three in- 
dividuals to 1C-E, two to 6C-E, and four to 
8C-E. A follow-up CVA (table 8) used these 
assignments to depict the pattern of morpho- 
logical variation and skin exchanges among 
the nine specimens. In addition to two color 
pattern characters disrupted by skin trans- 
plants, the PSC character did not meet the 
minimum tolerance criterion of the CVA 
model and was also omitted from the anal- 
ysis. Despite the reduced number of charac- 
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ters, the classification success was 100% 
within each of the three morphological sub- 
groups: 1C-E transplants, 6C-E transplants, 
and 8C-E transplants. Reciprocal skin trans- 
plants had been made (and accepted) among 
representatives of all three morphological 
subgroups (fig. 10A). 

Cordes and Walker (2003) also demon- 
strated unequivocal histocompatibility be- 
tween all lizards of pattern classes C-E and 
D tested. A CVA (table 9), using morpho- 
logical subgroups of C-E and representatives 
of pattern class D as a priori groups, revealed 
the pattern of morphological variation among 
the four groups, and substantiated that recip- 
rocal skin transplants had been made _ be- 
tween representatives of each morphological 
subgroup (1C-E, 6C-E, and 8C-E) and pat- 
tern class D (fig. 10B). 


PATTERNS OF MORPHOLOGICAL VARIATION 
AMONG PATTERN CLASSES C, D, AND E 
FROM SUMNER LAKE STATE PARK 


Different color pattern classes of A. tes- 
selata coexist at three localities—near the 
historic townsite of Higbee, Colorado, and in 
the vicinities of Conchas and Sumner Lakes, 
New Mexico. These assemblages provide 
valuable information on phenotypic diver- 
gence because shared environments increase 
the likelihood that their phenotypic differ- 
ences have a genetic basis. The Sumner Lake 
assemblage was particularly important be- 
cause this is the only locality known where 
pattern class E coexists with pattern classes 
C and D (Taylor et al., 1997). This assem- 
blage facilitated the resolution of more inclu- 
sive patterns of morphological variation in A. 
tesselata from eastern New Mexico. 

A CVA (table 10) was used to determine 
the morphological distinctiveness of the 
Sumner Lake pattern classes C (N = 45), D 
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Fig. 6. Pattern of multivariate morphological variation in three morphological subgroups of Aspi- 
doscelis tesselata of pattern class C-E (1C-E, N = 48; 6C-E, N = 18; and 8C-E, N = 13) from the 
vicinity of Conchas Lake State Park, San Miguel County, New Mexico. Samples of specimens of known 
genotype were used as a priori groups to make the original subgroup assignments; specimens used to 
make the assignments were omitted from the definitive CVA. A. Distributions derived from a CVA 
using nine meristic characters summarized in table 5. B. Distributions derived from a CVA of three 


color pattern characters summarized in table 6. 
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TABLE 6 
Discriminant Functions Used to Determine If 
Three Morphological Subgroups of 
Aspidoscelis tesselata of Pattern Class C—E 
from the Vicinity of Conchas Lake, San 
Miguel County, New Mexico, Could be 
Distinguished by Color Pattern Characters 


Charactera CV1 CV2 
DL-breaks — 1.288 —0.205 
PV-breaks 0.312 —0.254 
L-breaks 0.531 1.146 
Eigenvalue 0.114 0.021 


Total explained variation 84.5% 15.5% 


4Characters are defined in table 1. 


(N = 5), and E (N = 32). Individuals were 
assigned to a priori group by visual assess- 
ment of color pattern. The CVA gave an 
overall classification success of only 79%; 
however, except for 1 of 32 individuals of 


TABLE 7 
Discriminant Functions Used to Distinguish 
Three Morphological Subgroups of 
Aspidoscelis tesselata of Pattern Class C—E 
and Pattern Class D from the Vicinity of 
Conchas Lake, San Miguel County, New 














Mexico 
Character CV1 CV2 
GAB — 0.982 0.030 
LSG — 0.262 —0.048 
L-breaks 0.072 0.044 
GS 0.092 —0.056 
FP 0.143 —0.792 
COS 0.356 —0.518 
PSC 0.042 0.387 
DL-breaks 0.188 — 0.441 
PV-breaks —0.016 0.165 
Eigenvalue 3.643 0.911 
Total explained variation 76.1% 19.0% 





aCharacters are defined in table 1. 
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pattern class E misclassified to pattern class 
C, the other misclassifications were between 
pattern classes C and D (15 of 45 Sumner C 
assigned to Sumner D and | of 5 Sumner D 
assigned to Sumner C). The high number of 
misclassifications between C and D empha- 
sizes the meristic similarity of these two 
groups. 

The pattern of variation at Sumner Lake is 
a contrast between the meristic similarity of 
pattern classes C and D and the pronounced 
difference between these classes and Sumner 
E (fig. 11). Most individuals of Sumner E 
and Sumner C can be distinguished by visual 
inspection of the dorsal color pattern (e.g., 
Taylor et al., 1997: 864, fig. 1). However, a 
few individuals require subjective decisions 
regarding the degree of lateral barring (typ- 
ically less extensive in C, more extensive in 
E) and the pale vertebral line (typically pre- 
sent and relatively intact in C and absent or 
broken extensively in E). For example, cer- 
tain individuals of Sumner E have a relative- 
ly prominent vertebral line, thereby combin- 
ing a C-like feature with an E-like, disrupted 
lateral stripe (e.g., RU 9743 shown in fig. 
12A, B). The reverse combination also oc- 
curs; for example, RU 9507 (fig. 12D) was 
previously assigned to pattern class E be- 
cause of its disrupted vertebral field (Taylor 
et al., 1997, 2000). However, a CVA classi- 
fied this individual as pattern class C (P = 
0.90), with the assignment being influenced 
by the low number of stripe disruptions: 4 L- 
breaks, 0 DL-breaks, and 6 PV-breaks (fig. 
12E). In addition, certain individuals of Sum- 
ner C expressed unusual, partially divided 
vertebral lines, some just below the threshold 
of development seen in pattern class D from 
this locality (e.g., RU 9710 and 9725 shown 
in fig. 8B, C). As another example of con- 
founding color pattern elements, certain in- 
dividuals of Colorado D had partially divid- 
ed, but essentially single, vertebral lines 
(Taylor et al., 1996, fig. 1C, D: 256). 


-> 


Color pattern variation in Aspidoscelis tesselata of pattern class C-E from the vicinity of 


Conchas Lake State Park, San Miguel County, New Mexico. Morphological subgroup 6C-E: A (RU 
0002, 93 mm SVL); B (RU 0029, 96 mm SVL); morphological subgroup 1C-E: C (RU 0013, 95 mm 
SVL); D (RU 0030, 89 mm SVL); E (RU 0021, 95 mm SVL); morphological subgroup 8C-E: F (RU 


0027, 86 mm SVL). 
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COMPARISON OF PATTERN CLASSES C, C-E, 
E, AND E-C FROM EASTERN NEW MEXICO 


It appears likely that pattern classes C and 
C-E were derived from an ancestral pattern 
class, such as E or E-C. These are logical 
ancestral candidates because only popula- 
tions presently allocated to pattern class E 
are sympatric with the progenitor species of 
A. tesselata. Therefore, if E or E-C are older 
pattern classes, their lineages should have ac- 
crued more mutations than those of pattern 
classes C, C-E, and D. This possibility man- 
ifests itself in the color pattern variation 
found within and among groups traditionally 
assigned to pattern class E. A preview of this 
variation is provided by Walker et al. (1997a: 
240, fig. 4). Color pattern variation is ex- 
pected because clonal reproduction of selec- 
tively neutral, genetically modified pheno- 
types should include color pattern features. 
This would be evidenced by significant dif- 
ferences in color pattern characters between 
geographically proximate populations of the 
Same pattern class. As an example, the de- 
gree of lateral stripe disruption (L-breaks) 
was the principal character used by Zweifel 
(1965) to distinguish pattern classes C and E. 
Our sample of Luna E was not only signifi- 
cantly different in L-breaks from Sumner C 
(as expected), it also differed from Sumner 
E (P < 0.0005 for each comparison), al- 
though these samples came from sites sepa- 
rated by only 34 km. Paradoxically, there 
was no significant difference in L-breaks be- 
tween Luna E and each morphological sub- 
group of Conchas C-E (F313) = 2.225; P = 
0.09). 

We summarized color pattern variation by 
a CVA of color pattern characters (table 11), 
with samples of Conchas 1C-E (+D), 6C-E, 
and 8C-E, Sumner C (+D), Sumner E, Luna 
E, and Macho E-C serving as a priori groups. 
Because of insignificant multivariate differ- 
ences between Sumner D and Sumner C and 


Fig. 8. 
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between Conchas D and Conchas 1C-E, each 
of these pairs was pooled for this (and the 
following) analysis. The pattern of variation 
was a tight cluster of C-like individuals, a 
loose cluster of E-like individuals, and broad 
distributions of 1C-E, 6C-E, 8C-E, and Ma- 
cho E-C individuals across both C-like and 
E-like clusters (fig. 13A). 

We next used the seven a priori groups 
from the previous analysis in a CVA (table 
12) to depict a comprehensive (color pattern 
and scalation characters) pattern of morpho- 
logical variation. The inclusion of scalation 
characters enhanced the separation of Sum- 
ner E and Luna E from the other groups, with 
Conchas 6C-E and Macho E-C being distin- 
guished, as before, by the breadth of their 
morphological variation (fig. 13B). We sim- 
plified this two-dimensional pattern to one 
dimension by an additive tree phenogram 
(fig. 14) constructed from a symmetric ma- 
trix of Mahalanobis D? distances between 
sample centroids (table 13). Note that the re- 
semblance of Sumner D was with syntopic 
Sumner C (distance of 9.4) rather than with 
Conchas D (distance of 15.2), with smaller 
values reflecting greater similarity. The same 
pattern exists at Conchas Lake, where Con- 
chas D was most similar to syntopic Conchas 
1C-E (distance of 10.0) with which it shares 
the same GPI, EST2 genotype. Disregarding 
the Conchas D derivative, Conchas 1C-E and 
8C-E were morphologically most similar to 
Sumner C (fig. 14). 

Certain patterns of geographic variation 
suggest that a population resembling Macho 
E-C played an important role in the coloni- 
zation of northern habitats by A. fesselata. 
One example involves the three populations 
identified historically as pattern class E. De- 
spite the short geographic distance between 
the Sumner Lake and Puerto de Luna sites 
(fig. 1), Sumner E and Luna E resembled 
Macho E-C more closely than they resem- 


-> 


Color pattern variation in Aspidoscelis tesselata of pattern classes C and D from the vicinity 


of Sumner Lake State Park, De Baca County, New Mexico, and pattern class D from the vicinity of 
Conchas Lake State Park, San Miguel County, New Mexico. Sumner C: A (RU 9711, 97 mm SVL); B 
(RU 9710, 97 mm SVL); C (RU 9725, 93 mm SVL); Conchas D: D (RU 0004, 95 mm SVL); E (RU 
0043, 98 mm SVL); Sumner D: F (RU 9632, 97 mm SVL). Note the partial doubling of the vertebral 
stripe in individuals B and C (double vertebral lines being a characteristic of pattern class D). 
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bled each other (additive tree distances of 
17.0 for Sumner E—Macho E-C; 19.4 for 
Luna E—Macho E-C; 30.2 for Sumner E-— 
Luna E). In addition, Sumner C and the three 
morphological subgroups of Conchas C-E 
were also more similar to Macho E-C than 
to either of the geographically proximate 
populations of pattern class E. In fact, as the 
E-C designation suggests, certain individuals 
at Arroyo del Macho are C-like, with rela- 
tively intact pale stripes and vertebral lines 
(see Taylor et al., 2001: 16, fig. 4C, FE and 
compare these two specimens of pattern class 
E-C to representatives of pattern class C il- 
lustrated in Walker et al., 1997a: 239, fig. 3). 
Another example linking Macho E-C to 
northern groups was the strong morphologi- 
cal resemblance of Conchas 6C-E to Macho 
E-C rather than to syntopic clones of Con- 
chas 1C-E and 8C-E (fig. 14). 


STATUS OF PATTERN CLASS NEW MExIco D 


Each group of pattern class D, including 
Colorado D (Taylor et al., 1996), resembled 
a sympatric group of pattern class C or C-E 
rather than pattern class D from other local- 
ities. This suggests that pattern class D was 
either derived by mutation from its compan- 
ion population of pattern class C or C-E or 
that shared environments caused develop- 
mental convergence of meristic variation be- 
tween the two pattern classes at each site of 
sympatry. Evidence from the Sumner Lake 
pattern classes indicates that the environmen- 
tal alternative is unlikely. Despite shared 
habitats and activity periods, Sumner E is 
meristically distinct from Sumner C and 
Sumner D. Therefore, syntopy does not en- 
sure morphological similarity, and the close 
meristic resemblance between Sumner C and 
Sumner D is more likely based on genetic 
similarities between ancestral (C) and de- 
rived (D) clones. Because the sampling pro- 
cedure involved collecting individuals as 
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they were encountered, the ratio of 47 C:47 
E:5 D collected in 1995-1997 should ap- 
proximate the relative numbers of each pat- 
tern class at Sumner Lake. Although pattern 
class D is polyphyletic in the sense of having 
originated several times from different indi- 
viduals of pattern classes C and C-E, the pro- 
portion of its representatives at each locality 
does increase latitudinally from south to 
north: Sumner Lake (47 C:5 D = 10%), 
Conchas Lake (126 C-E:33 D = 21%), and 
Higbee (96 C:79 D = 45%, Taylor et al., 
1996). This trend could reflect different times 
of origin, different levels of success in dif- 
ferent habitats, or both. 

Zweifel (1965) discerned that pattern class 
C was distributed as a series of allopatric 
populations, whereas pattern class D was al- 
ways sympatric with pattern class C. Because 
of significant meristic differences between 
pattern class D from Higbee and Conchas 
Lake, Zweifel also hypothesized that pattern 
class D had originated twice from represen- 
tatives of pattern class C. Multivariate mor- 
phological evidence supports this hypothesis 
(Taylor et al., 1996). Based on identical ge- 
notypes for GPI, EST2, sACOH, and MPI 
and meristic similarities, we now identify 
Conchas clone IC-E as the likely progenitor 
of Conchas D and provide morphological ev- 
idence that Sumner D was similarly derived 
from Sumner C (fig. 14). Pattern class D is 
atavistic in the expression of broken, double 
vertebral lines—a color pattern feature of 
many individuals of A. tigris marmorata, the 
maternal progenitor species of A. tesselata. 


MODEL 1 FOR ORIGINS: MULTIPLE 
HYBRIDIZATIONS EXPLAIN THE 
VARIATION IN ASPIDOSCELIS 
TESSELATA IN EASTERN 
NEW MEXICO 


The only published evidence to explain the 
genetic and morphological variation in A. 
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Pattern of multivariate morphological variation in Aspidoscelis tesselata of subgroups 6C-E 


(N = 25), 1C-E (N = 58), 8C-E (N = 17), and pattern class D (N = 32) from Conchas Lake State 
Park, San Miguel County, New Mexico. Samples of specimens of known genotype were used as a priori 
groups to assign representatives of pattern class C-E to morphological subgroup. A. Distributions derived 
from a CVA using nine meristic characters (table 7). B. Centroids (means) of the distributions shown 


in part A. 
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TABLE 8 
Discriminant Functions Used in a CVA to 
Depict the Pattern of Morphological 
Variation Among Skin Transplant Specimens 
of Aspidoscelis tesselata of Pattern Class C—E 
from the Vicinity of Conchas Lake, San 
Miguel County, New Mexico 





Character CV1 CV2 





GAB —2.780 —0.514 
FP 2.255 —-0.990 
GS 2.208 — 0.404 
LSG 1.652 -0.591 
L-breaks 1.097 0.147 
COS —0.950 —0.510 
Eigenvalue 16.947 2.636 


Total explained variation 86.5% 13.5% 





4Characters are defined in table 1. 


tesselata from Conchas Lake is detailed in 
Parker (1979) and Parker et al. (1989). This 
hypothesis is based on allele combinations in 
A. tesselata being consistent with multiple 
hybridizations having taken place between its 
bisexual parental ancestors. Taken by itself, 
the allozyme diversity in the Conchas assem- 
blage could have been generated by multiple 
origins since there is no association between 
genotypes at the variable loci, and all four 
combinations of genotypes were detected in 
pattern class C-E (GPI —100/—96, EST2 
100/96; GPI —100/—96, EST2 96/96; GPI 
—99/—96, EST2. 100/96; *“GPI. —99/—96, 
EST2 96/96). For example, both the GPI 
—99 and —100 alleles appeared to be present 
in A. gularis septemvittata, and different hy- 
brid combinations could occur even within 
one clutch of eggs resulting from a single 
mating between representatives of A. tigris 
marmorata with a GPI —96/—96 genotype 
and A. gularis septemvittata having a GPI 
—100/—99 genotype. This would also be 
compatible with the morphological distinc- 
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tiveness of the clones marked by GPI —99/ 
—96 (Parker et al., 1989). The variation at 
EST2 is more problematical for this hypoth- 
esis. Because the 96-allele of EST2 was ap- 
parent only in A. gularis septemvittata and 
not in A. tigris marmorata, the origin of the 
EST2 96/96 genotype in the different GPI 
clones would require gene conversion (Hillis 
et al., 1991) or convergent mutation from the 
original F, hybrid genotype, EST2 100/96. 
However, if A. tigris marmorata previously 
had the 96-allele of EST2, or has it in low 
frequency or at localities not yet tested, then 
these genotypes could also have been pro- 
duced by different hybrid combinations of 
parental gametes. 

One problem with this model is biogeo- 
graphical—the necessity of a historical 
northward displacement of the current, geo- 
graphically restricted sympatry between the 
two progenitor species. Could the former 
range of A. gularis septemvittata have ex- 
tended farther north and overlapped that of 
A. tigris marmorata more extensively than it 
does today in southwestern Texas (Parker 
and Selander, 1976)? If so, hybridization 
could then have taken place closer to the pres- 
ent ranges of color pattern classes C, C-E, 
and D. Although the northernmost popula- 
tion of A. tigris marmorata in the Pecos Riv- 
er drainage (vicinity of Arroyo del Macho) 
is reasonably close to the Sumner Lake pat- 
tern classes, the northernmost population of 
A. gularis septemvittata is close to Deer 
Mountain, Hudspeth County, Texas (Dixon, 
2000). These populations are separated by 
approximately 400 linear km, and the Pecos 
River corridor is occupied by an intervening 
population of the more ecologically flexible 
A. gularis gularis (Degenhardt et al., 1996: 
217). A recent range shift of this magnitude 
is doubtful, and it is also unlikely to have 
taken place through the range of A. gularis 
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Fig. 10. Multivariate morphological variation among skin transplant specimens of Aspidoscelis tes- 
selata from Conchas Lake State Park, San Miguel County, New Mexico. Numbers are University of 
Arkansas, Department of Zoology (UADZ) accession numbers; lines identify reciprocal skin exchanges. 
Separate CVAs were run for each depiction (tables 8, 9). A. Transplants between individuals of three 
morphological subgroups of pattern class C-E. B. Transplants between representatives of three morpho- 
logical subgroups of pattern class C-E and New Mexico D. 
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TABLE 9 
Discriminant Functions Used in a CVA to 
Depict the Morphological Affinities of Skin 
Transplant Specimens of Aspidoscelis tesselata 
of Pattern Classes C—E and D from the 
Vicinity of Conchas Lake, San Miguel 
County, New Mexico 





Charactera CV1 CV2 
PSC — 1.394 —0.237 
LSG — 1.061 0.297 
COS 0.993 0.160 
L-breaks 0.763 0.541 
GAB —0.645 — 0.364 
FP — 0.265 0.990 
GS —0.071 0.427 
Eigenvalue 6.931 1.064 


Total explained variation 81.4% 12.5% 








Characters are defined in table 1. 


gularis without leaving evidence of intergra- 
dation. 


MODEL 2 FOR ORIGINS: 
POSTFORMATIONAL MUTATIONS 
EXPLAIN THE EVOLUTION OF 
ASPIDOSCELIS TESSELATA IN 
EASTERN NEW MEXICO 


One purpose of this study was to evaluate 
an alternative hypothesis to explain the 
source and pattern of genetic and morpho- 
logical variation in A. tesselata from eastern 
New Mexico. This hypothesis posits that the 
morphological diversity in A. tesselata is 
based on postformational genetic divergence 
(i.e., after establishment of parthenogenesis) 
following the origin of A. fesselata from a 
single F, hybrid zygote (Taylor et al., 1996; 
Walker et al., 1997b: 167, fig. 1; Taylor et 
al., 2000). 

This model is tied to three premises. The 
first assumes that histocompatibility among 
different individuals of Aspidoscelis tesselata 
reflects a single origin from one ancestral zy- 
gote. Although the genetic mechanism re- 
sponsible for the immune response in Aspi- 
doscelis has not been identified from exten- 
sive transplant experimentation (Cuellar and 
Smart, 1979), a rationale for interpreting his- 
tocompatibility as evidence of identical sets 
of histocompatibility alleles was provided by 
Cuellar (1976, 1977b), following Maslin 
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TABLE 10 
Discriminant Functions Used to Distinguish 
Aspidoscelis tesselata of Pattern Classes C, 
New Mexico D, and E from Sumner Lake 
State Park, De Baca County, New Mexico 


Characteré CV1 CV2 





L-breaks 0.897 0.299 
GAB 0.268 0.327 
PV-breaks — 0.238 ~0.671 
PSC —0.152 0.645 
COS —0.028 —0.483 
LSG —0.185 —0.288 
DL-breaks 0.191 0.094 
FP 0.279 —0.185 
GS 0.122 — 0.074 
Eigenvalue 7.203 0.078 


Total explained variation 98.9% 1.1% 





4Characters are defined in table 1. 


(1967). This rationale was used to explain 
histocompatibility among parthenogenetic 
individuals of A. cozumela and A. maslini 
(Hernandez-Gallegos et al., 1998). It was 
also used as evidence that two clones of A. 
laredoensis had been derived from separate 
zygotes (Abuhteba et al., 2000). Conversely, 
histoincompatibility evidence for separate 
hybridizations would be subject to error if 
mutations in histocompatibility alleles had 
occurred in clones derived from the same F, 
hybrid zygote. Possible examples of such in- 
dividuals are identified in A. laredoensis by 
Abuhteba et al. (2000). 

Although histocompatibility genetics in 
lizards is poorly understood, the demonstra- 
tion of histocompatibility for different color 
pattern classes of A. tesselata by Maslin 
(1967) and Cordes and Walker (2003) is a 
compelling argument for this model. If we 
assume that similar genetic systems are re- 
sponsible for histocompatibility in lizards 
and humans, then the genetic variation at the 
sACOH and MPI loci (table 3) is consistent 
with evidence that the clonal diversity at 
Conchas Lake State Park originated from a 
single F, hybrid zygote. There are two 
sACOH genotypes (be and ac) in A. tesselata 
at Conchas Lake. The origin of two of the 
three alleles can be accounted for because the 
c-allele is found in A. tigris marmorata, and 
the b-allele is found in A. gularis septemvit- 
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Pattern Class 


© Sumner C 
A Sumner E 
m Sumner D 


CV1 (98.9%) 


Fig. 11. Pattern of multivariate morphological variation among Aspidoscelis tesselata of pattern 
classes C (N = 44), E (NV = 32), and New Mexico D (N = 5) from the vicinity of Sumner Lake State 
Park, De Baca County, New Mexico. Canonical variate scores were derived from a canonical variate 


analysis using meristic characters identified in table 10. 


tata (the two progenitor species of A. tesse- 
lata, Neaves, 1969; Parker and Selander, 
1976; Dessauer and Cole, 1989; Dessauer et 
al., 1996). The common genotype be occurs 
in lizards of pattern classes C, D, E, and E- 
C, and clone VIC-E, and genotype ac is 
found only in certain individuals of clones 
IC-E and VIIIC-E. The a-allele has not been 
found in either bisexual progenitor taxon. 
The a-allele apparently does occur (with the 
b-allele) in A. gularis gularis (a close relative 
of A. g. septemvittata) and A. sexlineata 
(Dessauer and Cole, unpubl. data), so it is 
possible that the a-allele exists, undiscov- 
ered, in some extant population of A. g. sep- 
temvittata. 

The same argument can be made from the 
MPI alleles, with genotypes ab and ac found 


in A. tesselata from Conchas Lake. Again, 
two of the three alleles could have come 
from the progenitor taxa; the a-allele is found 
in A. tigris marmorata, the b-allele is found 
in A. gularis septemvittata, but the c-allele 
has not been found in either progenitor tax- 
on. The c-allele apparently does occur (with 
the b-allele) in A. gularis gularis (Dessauer 
and Cole, unpubl. data) so it too might also 
occur in some populations of A. g. septem- 
vittata. Therefore, if genotypic diversity at 
the sACOH and MPI loci were acquired from 
different alleles in separate fertilizations, it 
would require three separate F, hybrid zy- 
gotes from A. tigris marmorata (SACOH c-, 
MPI a-) X A. gularis (sACOH ab, MPI: bc) 
to generate the three sACOH, MPI clones at 
Conchas Lake (be ab; be ac; ac ac). Because 
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different combinations of histocompatibility 
alleles would be expected in each zygote 
(Maslin, 1967; Cuellar and Smart, 1979), 
skin graft rejection would then be expected 
between individuals within subgroups 1C-E 
(be ac and ac ac) and 8C-E (be ac and ac ac) 
and between individuals in these subgroups 
and those in 6C-E (be ab). A test of this hy- 
pothesis was provided by the reciprocal skin 
exchanges among all morphological sub- 
groups (fig. 10A). The complete lack of skin 
graft rejection is cause for rejecting the mul- 
tiple hybridizations hypothesis. 

The alternative hypothesis posits that this 
genotypic diversity resulted from point mu- 
tations at each of the two loci after the par- 
thenogenetic lineage had been established. 
The putative ancestral sACOH, MPI geno- 
type (bc ab) characterizes pattern class E-C 
from Arroyo del Macho and clone VIC-E 
from Conchas Lake. The alternative sACOH, 
MPI genotypes (be ac and ac ac), both in 
clones IC-E and VIIIC-E, could be accounted 
for by a single mutation at each locus. Sam- 
ples from other localities were not assessed 
for sACOH and MPI genotypes. 

The second premise assumes that evolu- 
tionary patterns in A. fesselata can be in- 
ferred from patterns of genotypic and mor- 
phological variation. Densmore et al. (1989) 
concluded, from the low level of sequence 
divergence in mtDNA, that A. tesselata had 
originated very recently. The same conclu- 
sion had been reached by Parker and Selan- 
der (1976) based on the relatively low ge- 
notypic variation found throughout the range 
of this species. In addition, A. tesselata lacks 
a mechanism (sexual reproduction) for inter- 
nal cohesion of variation (Dessauer and Cole, 
1986). Therefore, if A. tesselata originated 
from a single hybridization event (one F, zy- 
gote), spontaneous mutation would produce 
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TABLE 11 
Discriminant Functions, Derived from Color 
Pattern Characters, Used to Distinguish 
Aspidoscelis tesselata of Pattern Classes C—E 
(Three Morphological Subgroups) from 
Conchas Lake State Park, San Miguel 
County; Pattern Classes C and E from 
Sumner Lake State Park, De Baca County; 
Pattern Class E from Puerto de Luna, 
Guadalupe County; and Pattern Class E—C 
from Arroyo del Macho, Chaves County, All 
in New Mexico 
(Samples of pattern class D have been pooled 
with samples of putative local ancestral groups 
1C—E and Sumner C for this analysis.) 





Character CV1 CV2 
L-breaks 0.850 0.533 
DL-breaks 0.274 —~ 1.370 
PV-breaks — 0.066 0.879 
Eigenvalue 1.035 0.155 


Total explained variation 83.6% 12.5% 





4Characters are defined in table 1. 


a pattern of mosaic divergence into various 
clones. 

The third premise is based on the corre- 
spondence between meristic and genotypic 
variation identified in Parker (1979) and this 
study. We assume that genotypes of untested 
groups can be inferred from either geograph- 
ic proximity or morphological resemblance 
to groups of known genotype. For example, 
although genetic information was not avail- 
able for pattern class E from Sumner Lake, 
the GPI, EST2 genotypes are known for pat- 
tern class E from Fort Sumner, only 18 linear 
km southeast of the Sumner Lake sampling 
localities. All individuals from Fort Sumner 
(N = 11) belonged to genotypic clone VIE 
(GPI —99/—96 [Parker and Selander, 1976: 
798, table 5] and EST2 100/96 [Parker and 


Fig. 12. Color pattern variation in Aspidoscelis tesselata of pattern classes C and E from the vicinity 
of Sumner Lake State Park, De Baca County, New Mexico, and Puerto de Luna, Guadalupe County, 
New Mexico. Sumner E: A and B (RU 9743, 92 mm SVL) and C (RU 9748, 87 mm SVL); Sumner 
C: D and E (RU 9507, 89 mm SVL) (based on a visual assessment of the disrupted vertebral stripe, 
RU 9507 had been previously assigned to pattern class E. However, a CVA classified this individual as 
a pattern class C, with the relatively uninterrupted lateral stripe being important in this assignment.); 


Luna E: F (RU 0054, 95 mm SVL). 
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Selander, 1976: 797, table 4]). The Pecos 
River drainage provides a habitat corridor 
between Fort Sumner and Sumner Lake, and 
the general historical pattern of colonization 
of suitable habitats by A. tesselata has been 
from south to north (Parker and Selander, 
1976). In addition, Parker and Selander 
(1976) found GPI variation within individual 
pattern classes at only two widely separated 
localities—Conchas Lake, in the Canadian 
River drainage, and the south end of Ele- 
phant Butte Reservoir and Engle vicinity, 
both in the Rio Grande drainage. The geo- 
graphically restricted nature of genotypic 
variation in A. tesselata makes it likely that 
pattern class E at Sumner Lake State Park 
and Fort Sumner share the same GPI —99/ 
—96 genotype. This GPI genotype also char- 
acterizes clone VIC-E from Conchas Lake. 
Differences in the EST2 genotype have been 
found only in C-E from the Conchas Lake 
vicinity. 

We also lack genetic data for Sumner C, 
but we surmise that Sumner C and Sumner 
E have different GPI genotypes because of 
their pronounced morphological differences 
(table 2; figs. 12, 13). Therefore, the geno- 
type for Sumner C was inferred from its mor- 
phological resemblance to another group. 
Morphologically, Conchas 1C-E and 8C-E 
(both GPI —100/—96) most closely resemble 
Sumner C, and Sumner C most closely re- 
sembles Macho E-C (GPI —100/—96) (fig. 
14). This chain of resemblance suggests that 
Sumner C is also GPI —100/—96, which is 
the most common GPI genotype in this spe- 
cies (Parker and Selander, 1976). 


DIVERGENCE IN A. TESSELATA FROM EASTERN 
NEw MExIco 


BIOGEOGRAPHICAL PERSPECTIVE: The ances- 
tral pattern class in A. tesselata can be in- 
ferred from the geographic distributions of 
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TABLE 12 
Discriminant Functions, Derived from 
Scalation and Color Pattern Characters, Used 
to Distinguish Aspidoscelis tesselata Pattern 
Classes C—E (Three Morphological 
Subgroups) from Conchas Lake State Park, 
San Miguel County; Pattern Classes C and E 
from Sumner Lake State Park, De Baca 
County; Pattern Class E from Puerto de 
Luna, Guadalupe County; and Pattern Class 
E-C from Arroyo del Macho, Chaves County, 
all in New Mexico 
(Samples of pattern class D have been pooled 
with samples of putative local ancestral groups 
1C—E and Sumner C for this analysis.) 








Character@ CV1 CV2 
GAB 0.870 —0.513 
PV-breaks —0.186 0.031 
GS 0.116 —0.009 
COS —0.054 0.566 
FP 0.107 0.403 
L-breaks 0.277 0.512 
LSG —0.260 —0.379 
DL-breaks 0.249 0.266 
PSC —0.011 —0.123 
Eigenvalue 2.399 1.130 


Total explained variation 48.2% 22.7% 


4Characters are defined in table 1. 


its color pattern classes and progenitor spe- 
cies A. tigris marmorata and A. gularis sep- 
temvittata. Pattern class E (sensu Zweifel, 
1965) is the only pattern class that contacts 
the two progenitor species, with contact re- 
stricted to a few sites in northeastern Chi- 
huahua, Mexico, and southwestern Texas 
(Scudday, 1973; Scudday and Dixon, 1973; 
Parker and Selander, 1976). The present geo- 
graphic range of pattern class E (sensu Zwei- 
fel, 1965) extends northward from Trans-Pe- 
cos Texas, with the Rio Grande and Pecos 
River drainages being the principal dispersal 
routes (Cuellar, 1977a, fig. 3: 841). Pattern 


Patterns of multivariate variation among Aspidoscelis tesselata of morphological subgroups 


1C-E and Conchas D (N = 90), 6C-E (N = 25), and 8C-E (N = 17) from the vicinity of Conchas Lake, 
San Miguel County; pattern classes C and D (N = 50) and E (N = 32) from the vicinity of Sumner 
Lake State Park, De Baca County; pattern class E from Puerto de Luna, Guadalupe County (N = 16); 
and pattern class E-C from Arroyo del Macho, Chaves County (VN = 36), all from New Mexico. A. 
Canonical variate scores derived from a CVA of three color pattern characters identified in table 11. B. 
Canonical variate scores derived from a CVA using nine meristic characters identified in table 12. 
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TABLE 13 
Pairwise Mahalanobis Distances (D*) Among Centroids of Canonical Variate Scores for Samples 
of Aspidoscelis tesselata 6C—E, 1C—E, 8C—E, and Pattern Class D from the Vicinity of Conchas 
Lake State Park, San Miguel County; Pattern Classes C, D, and E from Sumner Lake State 
Park, De Baca County; Pattern Class E from Puerto de Luna, Guadalupe County; and Pattern 
Class E—C from Arroyo del Macho, Chaves County, all in New Mexico 
(This dissimilarity matrix was used to construct figure 14.) 



























6C-E IC-E 8C-E SumnerE Sumner C LunaE MachoE-C ConchasD Sumner D 
6C-E 0.00 
1C-E 16.00 0.00 
8C-E 12.92 8.47 0.00 
Sumner E 11.70 20.56 23.75 0.00 
Sumner C 5.47 10.22 5.95 19.40 0.00 
Luna E 16.54 19.95 19.07 19.49 16.86 0.00 
Macho E-C Tee 9.64 6.70 9.07 6.44 14.07 0.00 
Conchas D 13.02 2.11 5.02 24.07 7.58 19.24 9.55 0.00 
Sumner D 5.30 7.94 Vte 19.92 1.44 16.31 9.03 5.81 0.00 
te LUNA E classes C and D are first encountered near 
v the northern range boundary of pattern class 
E, at Sumner Lake State Park, De Baca 
County, New Mexico (Taylor et al., 1997). 
SUMNER_E Historically, the range of pattern class E ex- 
tended north to Santa Rosa Lake (named Los 
SUMNER_D Esteros Reservoir when constructed), Gua- 
dalupe County, New Mexico (documented 
SUMNER_C ; : 
=] by specimens in the Museum of Southwest- 
44 GONCHAS D ern Biology), but this population was evi- 







CE_1 
CE_8 


CE_6 


Fig. 14. Additive tree (phenogram), based on 
Mahalanobis D? distances (table 13), depicting 
meristic resemblance among nine groups of As- 
pidoscelis tesselata. Distances (similarities) be- 
tween groups are computed by adding lengths of 
nodes between groups of interest. Terminal nodes 
represent the nine groups, and internal nodes rep- 
resent horizontal distances between clusters. As 
an interpretation example, the resemblance be- 
tween Conchas 6C-E and Conchas 1C-E is 7.2 + 
3.2 + 2.2 + 2.8 + 5.6 = 21.0, while the resem- 
blance between Conchas 6C-E and Macho E-C is 
2 NOs lS ab bel" 13.2. 


dently extirpated when the reservoir basin 
filled (Brown, 1980; J. Applegarth, personal 
commun.). We have not seen representatives 
of A. tesselata at Santa Rosa Lake State Park 
during several recent visits to this locality, 
although J. M. Walker discovered in 1996 the 
population at Puerto de Luna, approximately 
13.5 km southeast of Santa Rosa. 

Scattered populations of pattern class C 
are found northeast of Sumner Lake in parts 
of northeastern New Mexico, northwestern 
Texas, southeastern Colorado, and south- 
western Oklahoma (Zweifel, 1965). There- 
fore, the relatively abrupt latitudinal shift 
from pattern class E to pattern class C may 
be related to mutations that preadapted pat- 
tern class C to the shorter growing seasons 
in the northern part of the range. The south- 
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ern range boundary of pattern class C (and 
D) at Sumner Lake is separated from sym- 
patric populations of the two bisexual pro- 
genitor species by approximately 520 linear 
km (Scudday, 1973; Scudday and Dixon, 
1973). Because pattern classes E-C and E are 
geographically interposed between the rang- 
es of the bisexual progenitor species and pat- 
tern classes C, C-E, and D, pattern classes 
E-C and E are ancestral pattern class candi- 
dates. The E-C alternative is particularly ap- 
pealing because of the breadth of its color 
pattern variation, although color pattern var- 
iants could also originate, similar to pattern 
class D, by mutation within established 
groups. 

MORPHOLOGICAL PERSPECTIVE: Most indi- 
viduals of A. tesselata can be assigned to col- 
or pattern class because of quasi-binary color 
pattern characters (Zweifel, 1965; Taylor et 
al., 1996; Walker et al., 1997a). Double, frag- 
mented vertebral lines distinguish most rep- 
resentatives of pattern class D from pattern 
classes C, C-E, E, and E-C, and the presence 
(Colorado D) or absence (New Mexico D) of 
supernumerary pale stripes distinguishes 
most individuals of the D-pattern classes 
(Taylor et al., 1996). Individuals of Colorado 
D lacking supernumerary stripes or having 
essentially intact vertebral lines are occasion- 
ally produced (Zweifel, 1965: 4, fig. 2D; 
Walker et al., 1997a: 238, fig. 2C). 

Because of the distinctive color pattern of 
New Mexico D, the assignment of individu- 
als to this pattern class is generally straight- 
forward. However, it is easy to demonstrate 
the reduced level of distinctiveness in Con- 
chas C-E and Macho E-C. Because the three 
morphological subgroups at Conchas Lake 
were most similar to either Sumner C (1C-E 
and 8C-E) or Macho E-C (6C-E), samples of 
Sumner C and Macho E-C were used as a 
priori groups in three discriminant analyses 
(DAs). Specimens of a different subgroup 
(1C-E, 6C-E, or 8C-E) were included as un- 
classified in each DA for classification to a 
priori group. The standardized discriminant 
functions (eigenvalue of 1.156) were 1.146 
for L-breaks, —0.607 for DL-breaks, and 
0.212 for PV-breaks. Classification success 
for the a priori groups was 96% for Sumner 
C (43 of 45 specimens assigned correctly) 
but only 73% for Macho E-C; 10 of 37 in- 


TAYLOR ET AL.: ASPIDOSCELIS TESSELATA 33 


dividuals were identified as belonging to 
Sumner C! With respect to the subgroups of 
pattern class C-E from Conchas Lake, the as- 
signments were (1) subgroup 1C: 36 individ- 
uals to Sumner C and 22 to Macho E-C, (2) 
subgroup 6C: 12 individuals to Sumner C 
and 13 to Macho E-C, and (3) subgroup 8C- 
E: 8 individuals to Sumner C and 9 to Macho 
E-C. This quantifies what is evident upon vi- 
sual inspection of Conchas Lake and Arroyo 
del Macho samples—use of pattern class C 
and E is misleading and oversimplifies color 
pattern variation in these populations. 

For example, the form at Arroyo del Ma- 
cho (Macho E-C) was referred to pattern 
class E by Zweifel (1965: 43), but differenc- 
es were noted by Taylor et al. (2001). The 
mix and range of color pattern expressions 
in Macho E-C were illustrated by a DA, us- 
ing the samples of Sumner E and Sumner C 
as a priori groups. These samples are appro- 
priate representatives of pattern classes C and 
E because of their high level of morpholog- 
ical distinctiveness (fig. 11). Individuals of 
Macho E-C were entered as unclassified for 
assignment to group. Because of the strong 
divergence between color patterns E and C 
at Sumner Lake, this model had a relatively 
high eigenvalue of 5.680, with standardized 
discriminant functions of 1.079 for L-breaks, 
0.091 for DL-breaks, and —0.254 for PV- 
breaks. There were few misclassifications in 
the a priori groups; 94% (30 of 32) of Sum- 
ner E and 100% (45 of 45) of Sumner C were 
classified correctly. However, 21 of 37 indi- 
viduals of Macho E-C were classified as 
Sumner C, with the remaining 16 classified 
as Sumner E. Although some populations of 
pattern class E lack evidence of a lateral 
stripe, Macho E-C has color pattern elements 
of both progenitor species—lateral barring 
from A. tigris marmorata and a lateral stripe 
from A. gularis septemvittata. The breadth of 
color pattern variation in Macho E-C is suf- 
ficient (fig. 13A) to have established the col- 
or pattern themes exhibited by the disjunct 
groups in northeastern New Mexico. 

GENETIC PERSPECTIVE: Ideally, a descrip- 
tion of the evolutionary history of A. tesse- 
lata would emerge from an initial under- 
standing of the processes generating the ge- 
netic and morphological variability in this 
clonally diverse species. Unfortunately, the 
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genetic bases of color pattern features and 
scalation differences are unknown, although 
there is evidence that elements of both are 
cloned (Dessauer and Cole, 1986; this study). 

Nevertheless, Densmore et al. (1989) pro- 
vide mtDNA evidence that pattern classes C, 
D, and E share one hybridization event in 
their past. This conclusion is consistent with 
the results of histocompatibility experiments 
(Maslin, 1967; Cordes and Walker, 2003). 
Therefore, if a single F, hybrid represented 
the origin of A. tesselata, its subsequent di- 
vergence into genetic and morphological 
clones should have a plausible explanation. 
As background information, genotype — 100/ 
—96 is the most common and widespread of 
the GPI alternatives (Parker and Selander, 
1976); therefore, it is the logical candidate 
for the ancestral combination of alleles. Ge- 
notype GPI —100/—96 characterizes popu- 
lations at Roswell (IE), Arroyo del Macho 
(IE-C), and Conchas Lake (IC-E and VIIIC- 
E). The presumed derivative, GPI —99/—96, 
is found at Fort Sumner (VIE) and Conchas 
Lake (VIC-E). 

The simplest fit of biogeographical, mor- 
phological, and genotypic data is provided 
by a model in which the origin and estab- 
lishment of derived color pattern classes and 
GPI clones followed a general south-to-north 
range expansion by colonists of GPI —100/ 
—96. Clone VIC-E from Conchas Lake 
shares the GPI —99/—96 genotype with VIE 
from Fort Sumner, yet the morphological af- 
filiations of both Conchas VIC-E and the 
Sumner Lake population of pattern class E 
are with clone IE-C from Arroyo del Macho. 
The mosaic nature of variation in different 
character systems (color pattern, scalation, 
proteins, karyotypes) is consistent with ran- 
dom mutations occurring after the establish- 
ment of parthenogenesis in A. tesselata. 

We discussed above two competing hy- 
potheses to explain the origin of variation in 
Aspidoscelis tesselata. Variation stemmed ei- 
ther (1) from different F, hybrid zygotes 
(separate fertilizations, whether involving the 
same individual parents and egg clutch or 
not), or (2) from mutations that occurred af- 
ter parthenogenesis was established. The hy- 
pothesis of multiple fertilization events is fa- 
vored by the fact that there are individuals 
of different color pattern classes (C, D, E) of 
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A. tesselata with vertebral dark field patterns 
that can be matched by color pattern varia- 
tion in A. tigris marmorata, and that certain 
variant alleles detected by protein electro- 
phoresis appear to exist as polymorphisms in 
the ancestral bisexual taxa. However, the ge- 
netic differences in these character states 
may be due to simple mutations or gene con- 
versions that can occur more than once, and 
they may be reversible. In addition, it is pos- 
sible that what appears on a gel to be the 
same allele product representing different liz- 
ards could actually be different allele prod- 
ucts having the same migration characteris- 
tics in the gel. 

We are impressed by the data produced by 
Maslin (1967) and Cordes and Walker (2003) 
demonstrating reciprocal histocompatibility 
among lizards of pattern classes C and E and 
among the genotypic clones of pattern class 
C-E and pattern class New Mexico D, re- 
spectively. If, as in humans, the histocom- 
patibility genes of Aspidoscelis involve mul- 
tiple loci with many independently assorting 
alleles (more than 50 alleles for one locus in 
humans; Gebhardt, 1996), this would suggest 
that A. tesselata originated from a single F, 
zygote. The rejection of skin transplants 
among 15 individuals from one population of 
bisexual A. figris is consistent with the as- 
sumption that lizards derived from different 
zygotes are histoincompatible (Cuellar and 
Smart, 1977). 


DEFICIENCIES OF THE PATTERN CLASS C AND 
E DESIGNATIONS 


The distinction between pattern classes C 
and E is essentially unambiguous at Sumner 
Lake, once color pattern ontogeny has been 
completed; for example, a CVA (table 11, 
fig. 13A) based on three color pattern char- 
acters classified all individuals correctly to 
two groups, C + D or E. Nevertheless, Zwei- 
fel (1965) provided examples that challenged 
the range-wide operationalism of pattern 
class E. This evidence included samples of 
pattern class E from disjunct sites in the Rio 
Grande drainage of New Mexico containing 
specimens with an intact lateral stripe (a 
characteristic of pattern class C). Zweifel 
(1965) also noted the unpredictable distri- 
bution of such variants in the geographic 
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range of pattern class E. Therefore, a dilem- 
ma exists for all who actually use color pat- 
tern criteria (rather than sampling locality) as 
the basis of color pattern assignments. Aban- 
doning pattern class C and E designations for 
divergent groups such as those at Sumner 
Lake is unnecessary. However, it is mislead- 
ing to refer to the non-D lizards at Conchas 
Lake as representing pattern class C or to the 
form at Arroyo del Macho as representing 
pattern class E. These populations comprise 
individuals with variable expressions of C 
and E color pattern elements and do not rep- 
resent syntopic assemblages of two distinct 
pattern classes. 


TAXONOMIC INTERPRETATION 


Small multivariate differences and shared 
habitats between Sumner C and Sumner D 
suggest that these pattern classes represent 
the same biological entity. This is also true 
for Conchas C-E and Conchas D, with the 
color patterns at each locality reflecting an- 
cestral (C and C-E) and derived (D) states. 
Therefore, New Mexico D is polyphyletic, 
having originated at least twice from differ- 
ent individuals in pattern classes C and C-E. 
Giving formal recognition to either group of 
New Mexico D is unwarranted. 

It is also apparent that the same species is 
represented by the three clones of pattern 
class C-E from the Conchas Lake vicinity. 
Although Conchas 1C-E, 6C-E, and 8C-E are 
genotypically and meristically distinctive, all 
three express the same color pattern varia- 
tion, an example of the mosaicism expected 
in clonal divergence. The three clones also 
appear to be ecologically equivalent based on 
sampling conducted (June 6—June 9, 2000) in 
one roadside habitat. The sample contained 
27 individuals exceeding 66 mm SVL, a size 
used in this study to exclude incompletely 
developed color pattern characters from the 
Statistical analyses. Of the 27 specimens, 11 
had been classified to 1C-E, 10 to 6C-E, and 
6 to 8C-E. Sampling variation can account 
for these numerical differences (Chi-square 
= 1.556; P = 0.46), indicating an absence of 
ecological segregation in this mesquite-dom- 
inated habitat. 

Among the groups studied, the most com- 
pelling case for postformational speciation 
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involves pattern classes C and E at Sumner 
Lake State Park. Besides color pattern dif- 
ferences, they are significantly different in 
meristic characters, mean body size (SVL) of 
reproductively mature individuals, and life- 
history characteristics. Sumner C and Sum- 
ner E occur at equivalent densities although 
Sumner C is larger than Sumner E and pro- 
duces larger clutches. Compensation is ap- 
parently achieved by certain individuals of 
Sumner E beginning reproduction in their 
second year, with Sumner C delaying repro- 
duction to year three (Taylor et al., 1997, 
2000). Sumner Lake State Park is of consid- 
erable interest as a natural experiment be- 
cause it represents an area where pattern 
classes E and C reach their northern and 
southern range limits, respectively. Do Sum- 
ner C and Sumner E reflect a postformational 
Speciation event, with independent trajecto- 
ries established by the separation of a toko- 
genetic clone vector (Wiley and Mayden, 
2000) or tokogenetic array (Frost and Hillis, 
1990) into two such entities? If so, a taxo- 
nomic decision regarding pattern class C is 
moot. The name A. fesselata is associated, 
through the selection of a neotype by Walker 
et al. (1997a), with Colorado D from Higbee. 
Therefore, unless Colorado D can be shown 
to represent a different biological entity 
(Walker et al., 1997b), pattern class C will 
continue to share the name A. tesselata with 
its Colorado D derivative. 

Groups such as Conchas C-E and Macho 
E-C demonstrate the discordance of color 
pattern features and geographic range bound- 
aries between entities traditionally regarded 
as pattern classes C and E. These two groups 
exemplify the challenges that mosaic patterns 
of genetic and phenotypic variation impose 
on the range-wide utility of color pattern 
classes and on the formal taxonomic recog- 
nition of pattern class “‘E”’ and its variants. 
The alternative to taxonomic restructuring is 
not without merit—a perspective of A. fes- 
selata as a single species, conceptually im- 
portant for demonstrating the effectiveness of 
a generalized genotype (Taylor et al., 2001), 
and subsequent random mutations, to meet 
evolutionary challenges in a parthenogenetic 
entity of recent origin. 
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APPENDIX 1 


SPECIMENS EXAMINED 


Specimens are referred to by their individual 
catalog numbers and initials of their repositories 
as follows: AMNH (American Museum of Natu- 
ral History), UADZ (University of Arkansas De- 
partment of Zoology [EDP = specimens collected 
by E.D. Parker, Jr; GM = UADZ specimens col- 
lected by Glenn J. Manning]); RU (Regis Univer- 
sity Department of Biology). Individuals classi- 
fied by CVA to morphological subgroup of pat- 


tern class C-E are identified by designations 1C- 
E, 6C-E, and 8C-E and individuals karyotyped for 
this paper with (k). Specimens lacking subgroup 
designations did not meet the minimum SVL cri- 
terion for inclusion in the CVA used to assign 
individuals to subgroup. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS C-E) 
WITH ELECTROPHORETIC DATA: New Mexico, San 
Miguel County, around picnic tables south of 
Conchas Dam: clone IC-E: UADZ 5418, 5419, 
5421-5423 (= EDP 836, 837, 839-841); clone 
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VIC-E: UADZ 5415-5417, 5426, 5427 (= EDP 
830, 831, 834, 845, 849); clone VIIC-E: UADZ 
5420, 5424, 5425, 5428 (= EDP 838, 842 [842 
was identified as a pattern class D in Parker and 
Selander, 1976; Parker, 1979; and Parker et al., 
1989], 844, 850. Conchas Lake at South State 
Park campground: clone VIC-E and sACOH (bc), 
MPI (ab): AMNH R-136875, R-136878 (k); clone 
IC-E and sACOH (bc), MPI (ac): AMNH R- 
123029, R-136877; clone IC-E and sACOH (ac), 
MPI (ac): AMNH R-123033 (k), R-136876, R- 
136879. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS C-E), 
LAB HATCHLINGS WITH ELECTROPHORETIC DATA: 
AMNH R-123026 (from egg produced by R- 
123025), R-123035 (from egg produced by R- 
123033); 

ASPIDOSCELIS TESSELATA (PATTERN CLASS C-BE), 
LAB HATCHLINGS LACKING ELECTROPHORETIC 
Data: AMNH R-123027, R-123028 (from eggs 
produced by R-123025), R-123030—123032 (from 
eggs produced by R-123029), R-123034, R- 
123036 (from eggs produced by R-123033). 

ASPIDOSCELIS TESSELATA (PATTERN CLASS C-BE), 
SKIN TRANSPLANT SPECIMENS: New Mexico, San 
Miguel County, 3 km south of Conchas Dam or 
east of south end of levee, along east side Hwy 
129 (Hwy 433): morphological subgroup 6C-E: 
UADZ 3243, 3449; morphological subgroup 1C- 
E: UADZ 3237, 3241, 3444; morphological sub- 
group 8C-E: UADZ 3240, 3246, 3247, 3448. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS C-—E) 
LACKING ELECTROPHORETIC DATA: New Mexico, 
San Miguel County. Along Army Corps of En- 
gineer road east of and paralleling Hwy 433 (the 
entrance road to Conchas Lake State Park), ap- 
proximately 1 km from the junction of this road 
and Hwy 433 (east of the entrance to Central Rec- 
reation Area): RU 0001 (8C-E), 0002 (6C-E), 
0003 (6C-E), 0008 (8C-E), 0009 (6C-E), 0010 
(8C-E), 0011 (1C-E), 0012 (6C-E), 0013 (1C-E), 
0014 (1C-E), 0015 (6C-E), 0016, 0017 (8C-E), 
0018 (1C-E), 0019 (6C-E), 0020 (1C-E), 0021 
(1C-E), 0022 (1C-E), 0027 (8C-E), 0028 (1C-E), 
0029 (6C-E), 0030 (1C-E), 0031 (6C-E), 0032 
(6C-E), 0033 (8C-E), 0037, 0038 (6C-E), 0039 
(1C-E), 0040 (6C-E), 0041 (1C-E). North side Ca- 
nadian River, 200-300 m below dam: GM 110 
(1C-E), 111 (1C-E), 112 (1C-E); south side Ca- 
nadian River, near one-way road into park, below 
dam: GM 113, 114; south side Canadian River, 
plateau above river: GM 116 (8C-E); south side 
Canadian River, parking lot below dam: GM 119, 
120 (1C-E); south side of Canadian River, South 
Recreation Area: GM 129; Central Recreation 
Area fronting lake: GM 136 (6C-E), 137 (1C-E), 
138, 139; north side Canadian River, north of Park 
Office and North Recreation Area near Boy Scout 
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camp: GM 142, 143 (1C-E), 145 (1C-E), 146 (1C- 
E); north side of Canadian River, near (below) 
North Recreation Area: GM 147; north side of 
Canadian River, near (above) North Recreation 
Area: GM 152 (1C-E). Conchas Lake at South 
State Park campground: EDP UADZ 843 (1C-B), 
844 (1C-E), 848; UADZ 3214, 3215 (1C-B), 
3217, 3218 (1C-E), 3219 (6C-E), 3220 (6C-E), 
3221, 3222 (1C-E), 3223 (8C-E), 3224 (1C-E), 
3226 1C-E), 3227 (6C-E), 3229, 3231 (6C-E), 
3232 (1C-E), 3233 (1C-E), 3234 (1C-E), 3236, 
3239 (1C-E), 3245 (8C-E), 3432 (8C-E), 3433, 
3434 (1C-E), 3435 (6C-E), 3437, 3438 (1C-E), 
3439 (1C-E), 3442, 3443 (1C-E), 3744 (1C-B), 
3745 (1C-E), 3746 (1C-E), 3747 (1C-E), 3748 
(1C-E), 3749 (1C-E), 4158 (1C-E), 4159 (1C-E), 
4160, 4161, 4162 (1C-B), 4163-4167, 4168 (1C- 
E), 4169 (1C-E), 4170 (1C-E). South State Park 
campground: AMNH_ R-119545 (8C-E), R- 
123025 (6C-E; k), R-123050 (6C-E), R-123051 
(8C-E; k). 

ASPIDOSCELIS TESSELATA (PATTERN CLASS C) 
LACKING ELECTROPHORETIC DATA: New Mexico, 
De Baca County, Sumner Lake State Park, West- 
side Campground area: RU 9507, 9651; terrace on 
west side Pecos River via road south of Hwy 203 
to West River Picnic area, then on flats north of 
unimproved road providing access to area to the 
west: RU 9735; east side of Pecos River and Sum- 
ner Lake from arroyo on south side of Hwy 203 
(southwest of the road to Eastside Campground), 
then 5.0 km north of Hwy 203 along road to East- 
side Campground and continuing east and north 
to primitive area: RU 9518, 9521, 9533-9537, 
9539-9541, 9602-9604, 9608-9611, 9614, 9615, 
9626, 9627, 9629, 9630, 9633, 9637, 9638, 9640— 
9642, 9647, 9709-9712, 9725-9729, 9737-9740, 
9745, 9746. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS D) 
WITH ELECTROPHORETIC DATA: New Mexico, San 
Miguel County, Conchas Lake State Park, around 
picnic tables south of Conchas Dam: clone ID: 
UADZ 5404-5414 (= EDP 826-829, 832, 833, 
835 [835 was identified as a pattern class C in 
Parker and Selander, 1976 and Parker, 1979], 843, 
846-848). Conchas Lake at South State Park 
campground: clone sACOH (bc), MPI (ac): 
AMNH R-123042, R-136880. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS D), 
LAB HATCHLINGS WITH ELECTROPHORETIC DATA: 
AMNH R-123038 (from egg produced by R- 
123037) R-123043, R-123044 (from eggs pro- 
duced by R-123042). 

ASPIDOSCELIS TESSELATA (PATTERN CLASS D), 
LAB HATCHLINGS LACKING ELECTROPHORETIC 
Data: AMNH R-123039-R123041 (from eggs 
produced by R-123037); AMNH R-123045, R- 
123046 (from eggs produced by R-123042). 
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ASPIDOSCELIS TESSELATA (PATTERN CLASS D), 
SKIN TRANSPLANT SPECIMENS: New Mexico, San 
Miguel County, 3 km south of Conchas Dam or 
east of south end of levee, along east side of Hwy 
129 (Hwy 433): UADZ 3238, 3242, 3244, 3450. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS D) 
LACKING ELECTROPHORETIC DATA: New Mexico, 
San Miguel County. Along Army Corps of En- 
gineer road east of and paralleling Hwy 433 (the 
entrance road to Conchas Lake State Park), ap- 
proximately 1 km from the junction of this road 
and Hwy 433 (east of the entrance to Central Rec- 
reation Area): RU 0004—0007, 0023-0026, 0034— 
0036, 0042—0046. South side of Canadian River, 
200-300 m below spillway, very close to river: 
GM 121, 122, 141. South of Canadian River Cen- 
tral Recreation Area of park, fronting lake: GM 
134. South State Park campground: AMNH R- 
114236 (k), R-114237 (k), R-119546, R-123037 
(k), R-123042, R-123047-123049, R-123052. 
New Mexico, De Baca County, Sumner Lake 
State Park, east side of Pecos River and Sumner 
Lake from arroyo on south side of Hwy 203 
(southwest of the road to Eastside Campground), 
then 5.0 km north of Hwy 203 along road to East- 
side Campground and continuing east and north 
to primitive area: RU 9517, 9612, 9632, 9741; 
terrace on west side Pecos River via road south 
of Hwy 203 to West River Picnic area, then on 
flats north of unimproved road providing access 
to area to the west: UADZ 5744. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS E) WITH 
ELECTROPHORETIC DATA: New Mexico, De Baca 
County, Fort Sumner city dump: EDP 309-319. 
New Mexico, Chaves County, Roswell city dump: 
EDP 631-640. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS E) 
LACKING ELECTROPHORETIC DATA: New Mexico, 
De Baca County. Sumner Lake State Park, West- 
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side Campground area: RU 9508-9516, 9519, 
9520, 9616; terrace on west side Pecos River via 
road south of Hwy 203 to West River Picnic area, 
then on flats north of unimproved road providing 
access to area to the west: RU 9605, 9736; east 
side of Pecos River and Sumner Lake from arroyo 
on south side of Hwy 203 (southwest of the road 
to Eastside Campground), then 5.0 km north of 
Hwy 203 along road to Eastside Campground and 
continuing east and north to primitive area: RU 
9522, 9523, 9538, 9542, 9543, 9601, 9605, 9613, 
9625, 9628, 9631, 9634, 9635, 9639, 9643-9646, 
9648, 9713-9716, 9730-9732, 9742-9744, 9747— 
9750. Fort Sumner, east end of Hwy 60-84 by- 
pass at base of elevated railroad grade and in ad- 
jacent mesquite assemblage: UADZ 5688-5690, 
5699. Fort Sumner, at Hwy 60—84 bypass, in mes- 
quite at entrance to city dump (Parker and Selan- 
der [1976] collecting site): UADZ 6308. New 
Mexico, Guadalupe County, south along sandy 
road, southeastern edge of Puerto de Luna: RU 
0048-0073. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS E-C) 
WITH ELECTROPHORETIC DATA: New Mexico, 
Chaves County, north side Arroyo del Macho, 22 
km north on Hwy 285 from junction with Hwy 
70 north of Roswell, New Mexico, then 0.8 km 
east on Eden Valley Road: AMNH R-145142— 
145144, 146629, 146630 (laboratory hatchling; 
offspring of either AMNH R-146629 or 146631), 
14663 1—-146633, 146636, 146638, 146639. 

ASPIDOSCELIS TESSELATA (PATTERN CLASS E-C) 
LACKING ELECTROPHORETIC DATA: New Mexico, 
Chaves County, north side Arroyo del Macho, 22 
km north on Hwy 285 from junction with Hwy 
70 north of Roswell, New Mexico, then 0.8 km 
east on Eden Valley Road: AMNH R-146612— 
146628, R-146634, R-146635, R-146637, R- 
146640-146649, R-147547. 
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